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Abstract 

Edwards' functional integral approach to the statistical mechanics of polymer liq- 
uids is amenable to a diagrammatic analysis in which free energies and correlation 
functions are expanded as infinite sums of Feynman diagrams. This analysis is shown 
to lead naturally to a perturbative cluster expansion that is closely related to the 
Mayer cluster expansion developed for molecular liquids by Chandler and cowork- 
ers. Expansion of the functional integral representation of the grand-canonical par- 
tition function yields a perturbation theory in which all quantities of interest are 
expressed as functionals of a monomer-monomer pair potential, as functionals of 
intramolecular correlation functions of non- interacting molecules, and as functions 
of molecular activities. In different variants of the theory, the pair potential may 
be either a bare or a screened potential. A series of topological reductions yields a 
renormalized diagrammatic expansion in which collective correlation functions are 
instead expressed diagrammatically as functionals of the true single-molecule cor- 
relation functions in the interacting fluid, and as functions of molecular number 
density. Similar renormalized expansions are also obtained for a collective Ornstein- 
Zernicke direct correlation function, and for intramolecular correlation functions. 
A concise discussion is given of the corresponding Mayer cluster expansion, and 
of the relationship between the Mayer and perturbative cluster expansions for liq- 
uids of flexible molecules. The application of the perturbative cluster expansion to 
coarse-grained models of dense multi-component polymer liquids is discussed, and a 
justification is given for the use of a loop expansion. As an example, the formalism 
is used to derive a new expression for the wave-number dependent direct correlation 
function and recover known expressions for the intramolecular two-point correlation 
function to first order in a renormalized loop expansion for coarse-grained models 
of binary homopolymer blends and diblock copolymer melts. 
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1 Introduction 



Diagrammatic expansions have played an important role in the development 
of the theory of classical fluids. They have proved useful both as a tool for the 
development of systematic expansions in particular limits, and as a language 
for the analysis of proposed approximation schemes. Analysis of Mayer cluster 
expansions of properties of simple atomic liquids was brought to a very high 
degree of sophistication by the mid-1960s [1,2,3]. Corresponding Mayer clus- 
ter expansions were later developed for interaction site models of both rigid 
[4] and non-rigid molecules [5,6,7] by Chandler and coworkers. Throughout, 
the development of both atomic and molecular liquid state theory has been 
characterized by an interplay between diagrammatic analysis [8,9] and the 
development of integral equation approximations [10,11,12,13,14,15,16]. 

A field-theoretic approach that was introduced by Edwards [17,18] has been 
used in many studies of coarse-grained models of polymer Uquids. Edwards' ap- 
proach relies upon an exact transformation of the partition function of a poly- 
mer fluid from an integral with respect to monomer positions to a functional 
integral of a fluctuating chemical potential field. A saddle-point approximation 
to this integral leads [19] to a simple form of mean field theory, which further 
reduces to a form of Flory-Huggins theory in the case of a homogenous mixture. 
The first applications of this approach were studies of excluded volume effects 
polymer solutions, by Edwards, Freed, and Muthukumar [17,18,19,20,21,22]. 
Several authors have since used Edwards' formalism to study corrections to 
mean field theory in binary homopolymer blends and in block copolymer melts 
[24,25,26,27,28,29,30,31,32,33,34,35]. All of these studies of dense liquids have 
been based upon some form of Gaussian, or, equivalently [23], one-loop ap- 
proximation for fiuctuations about the mean field saddle-point. This formalism 
is also the basis of a numerical simulation method developed by Fredrickson 
and coworkers [36,37], in which the functional integral representation of the 
partition function is sampled stochastically. 

Edwards' approach lends itself to the use of standard methods of perturbative 
field theory, including the use of Feynman diagrams. By analogy to experience 
with both Mayer cluster expansions in the theory of simple liquids, and of 
Feynman diagram expansions in statistical and quantum field theories, one 
might expect it to be possible to develop systematic rules for the expansion of 
the free energy and various correlation functions as well-defined infinite sets 
of Feynman diagrams. Such rules have, however, never been developed for 
Edwards' field-theoretic approach with a level of generality or rigor comparable 
to that attained long ago for either Mayer clusters or statistical field theory. 
This paper attempts to rectify this, while also exploring connections between 
different diagrammatic approaches to liquid state theory. 
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The analysis given here starts from a rather generic interaction site model 
of a flTiid of non-rigid molecules. Molecules are comprised of point-like parti- 
cles (referred to here as "monomers" ) that interact via a pairwise additive 
two-body interaction, and via an unspecified intramolecular potential among 
monomers within each molecule. The diagrammatic expansion that is obtained 
by applying the machinery of perturbative field theory to a functional inte- 
gral representation of the grand partition function is a form cluster expansion, 
which is referred to here as an interaction site perturbative cluster expansion. 
Terms in this expansion are conveniently represented in terms of diagrams of 
bonds and vertices, in which vertices represent multi-point intramolecular cor- 
relation functions. In different variants of the theory, the bonds may represent 
either the screened interaction identified by Edwards [17,18] or the bare pair 
potential. 

The expansion in terms of the bare potential is shown to be a particularly close 
relative of the interaction site Mayer cluster expansion developed for fluids of 
non-rigid molecules by Chandler and Pratt [5]. The main differences are dif- 
ferences in the diagrammatic rules that arise directly from the association of 
a factor of the pair potential rather than the corresponding Mayer /-function 
with each bond. A self-contained derivation of the Chandler-Pratt Mayer clus- 
ter expansion for a molecular liquid is given in Sec. 22. The derivation of the 
Mayer cluster expansion given here follows a line of reasoning closely analogous 
to that used in the theory of simple liquids, which starts from an expansion 
of the grand-canonical partition function. The derivation is somewhat more 
general and (arguably) more direct than that given by Chandler and Pratt. 

One result of the present analysis that does not seem to have an analog in the 
Mayer cluster analysis of Chandler and Pratt is an expansion of a generalized 
Ornstein-Zernicke direct correlation function for a fluid of flexible molecules. 
This is presented in Sec. 17. 

Throughout, the analyses of both perturbative and Mayer cluster expansions 
proceed by reasoning that is, as much as possible, analogous to that given 
for atomic liquids by Morita and Hiroike [1], Stell [2], and Hansen and Mac- 
Donald [3]. Diagrammatic rules for the calculation of correlation functions 
are derived by functional differentiation of an expansion of a grand canon- 
ical thermodynamic potential with respect to fields conjugate to monomer 
concentrations. Several renormalized expansions are obtained by topologically 
reductions roughly analogous to those applied previously to atomic fiuids. 
Because the diagrams used to construct perturbative cluster expansions for 
fiuids of non-rigid molecules are different than those used in either the Mayer 
cluster expansion for atomic liquids or in other applications of statistical field 
theory, the main text is supplemented by a discussion in appendices C and 
D of the symmetry numbers needed to calculate combinatorical prefactors for 
such diagrams, and in appendices E and F by several lemmas about diagrams 
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that are needed to justify various topological reductions. The required lemmas 
are all generalizations of those given by Morita and Hiroike for fluids of point 
particles. 

Once the close relationship between the perturbative and Mayer cluster ex- 
pansions is appreciated, it reasonably to ask for what classes of problems a 
perturbative expansion might be useful. A perturbative description of a clas- 
sical fluid is useful only for models with relatively soft or long-range pair 
interactions. It is thus clearly not a useful starting point for treating the harsh 
repulsive interactions encountered in any atomistic model of a dense liquid. As 
discussed in Sec. 19, a perturbative expansion is, however, potentially useful 
for the study of very coarse-grained models of dense multi-component polymer 
liquids [35] . In such models, in which each monomer is a soft "blob" that repre- 
sents a subchain of many chemical monomers, the effective interaction between 
blobs is much softer and of much longer range than the interactions between 
atoms or chemical monomers. It is shown here that a loop expansion of the 
perturbative diagrammatic expansion for such a model yields an asymptotic 
expansion about the mean fleld theory when the coarse-grained monomers are 
taken to be large enough so as to strongly overlap. The small parameter in this 
expansion is the ratio of the packing length of the melt (which is independent 
of the degree of coarse-graining) to the size of a coarse-grained monomer. 

As an example, in Sec. 20 the formalism is used to derive expressions for 
the two-point intramolecular correlation function and the direct correlation 
function in both binary homopolymer mixtures and block copolymer melts 
to first order in a renormalized loop expansion. The resulting expressions are 
compared to those obtained by related methods in several previous studies. 

Notwithstanding its title, very little in this article is specific to high molecular 
weight polymers: All of the results, except those of Sec. 19, are formally appli- 
cable to any interaction site model of a classical fiuid of non-rigid molecules. 
The analysis is nonetheless reasonably described as a theory of polymer liq- 
uids both for sociological reasons, because Edwards' approach has been used 
primarily by polymer physicists for the study of coarse-grained models of poly- 
mer liquids, and for practical reasons, because this approach seems to be best 
suited for this purpose. 



2 Model and Definitions 

We consider a mixture of S molecular species, each of which is constructed 
from a palette of C types of monomer. Each molecule of species a, where 
a = 1, . . . , S, contains Nai monomers of type i, where i = 1, . . . ,C. The total 
number of molecules of type a is Ma in a system of volume V, giving a number 
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density pa = {Mg) /V . The following conventions are used throughout for the 
names and ranges of indices: 



a,b = 


1,.. 


..,S 


molecule species 


hi = 


1,.. 


..,c 


monomer types 


m — 


1,.. 




molecules 


s — 


1,.. 


N ■ 


monomer sites 



(1) 



Unless otherwise stated, summations over indices are taken over the ranges 
indicated above, and summation over repeated indices is implicit. 

Let Tiami{s) be the position of monomer s of type i on molecule m of species 
a. We define monomer density fields 

Cami {r) = ^S(r- Rami («) ) 

s 

Ci{r) = J2cami{r) . (2) 

am 

These give the number concentrations for monomers of type i on a specified 
molecule m of type a, and for all monomers of type i, respectively. 



2.1 Model 



We consider a class of models in which the total potential energy of a molecular 
liquid is given by a sum 

U^Uo + Uint + Uea^t ■ (3) 



Here, Uo is an intramolecular potential energy. This might be taken to be a 
Gaussian stretching energy in a coarse-grained model of a polymer chain, but 
its exact form may be left unspecified. The energy Uint is a pairwise-additive 
interaction between monomers, of the form 

Uint ^lT.JdrJdr' U,j{r - r%{r)cj{r') (4) 



in which Uij{r — r') is the interaction potential for monomers of types i and 
j. The external potential is of the form 

Uext = - X] j^^ hi{^)ci{r) , (5) 
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where hi{r) is an external field conjugate to Ci(r). The external field h is intro- 
duced solely as a mathematical convenience, so that expressions for correlation 
functions may be derived by functional differentiation of the grand potential 
with respect to h. 

The canonical partition function for a system with a specified number of 
molecules of each type a, as a functional of the multi-component field h, is 
denoted 

1 _„ 



where M denotes a set of values M = {Mi, . . . , Mc} for all species, and D[R] 
denotes an integral over the positions of all monomers in the system. Here 
and hereafter, we use units in which /c^T = 1. The grand-canonical partition 
function is 

M 



where is a chemical potential for molecular species a, which corresponds to 
an activity Aa = e'^", and J2m denotes a sum over Ma > for all a. The second 
equality in Eq. (7) introduces the notation Tr[- ■ ■] to represent integration over 
all distinguishable sets of monomer positions and summation over M^. 



2.2 Notational Conventions 



Functions of n monomer positions and n corresponding monomer type indices 
may be expressed either using a notation 

i^i(l) = i^n...i„(ri,...,r„) (8) 



in which i = {ii, i2, ■ ■ ■ , in} and r = {ri, r2, . . . , r„} denote lists of monomer 
type and positions arguments, respectively, or in an alternative notation 



F(l,...,n) = F,,...,„(ri,...,r„) (9) 

in which an integer argument 'j' represents both a position rj and a monomer 
type index ij. 
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The notation H — G * F, when apphed to two-point functions, represents a 
convolution 



in which / d{j) = J2j Jdrj denotes both integration over a coordinate and 
summation over the corresponding type index. Similarly, h*cis shorthand for 
an integral / d{l)h{l)c{l). The function F^^{r,r') denotes the integral oper- 
ator inverse of a two-point function Fij{r, r'), which is defined by an integral 
equation by F * F~^ — S where 5(1,3) = 5i^igS{ri — r^). 

Coordinate-space and Fourier representations of functions of several variables 
may be used essentially interchangeably in most formal relationships. Fourier 
transforms of fields such as the monomer density Cj(r) are defined with the 
convention Cj(k) = Jdr e~*^ '"Cj(r), while transforms of functions of two or 
more monomer positions and type indices, such as Uij {ri, r2), are defined with 
the convention 



where / d"'r = J dri ■ ■ ■ J rir„, and where k = {ki, k2, . . . , k„} denotes a list of 
n wavevector arguments. The prefactor of 1/V in Eq. (11) guarantees that the 
transform Fi(k) of a translationally invariant function will approach a finite 
limit as y — > oo if ki -|- k2 -|- • • • + k„ = (and must vanish otherwise). When 
this convention matters, an n-point function -Fi(k) defined by convention (11) 
will be referred to as a normalized Fourier transform of i^i(r), and V^Fj(k) as 
the unnormalized transform. The normalized transform of a translationally 
invariant two-point function, such as Uij{r — r'), may be expressed as a func- 
tion f/ij(k) = [/y(k, — k) of only one wavevector. The normalized transform 
FjJ^(k) = Fjj^(k, k') of the inverse of a translationally invariant function F is 
simply the matrix inverse of Fjj(k). 

2.3 Collective Correlation and Cluster Functions 

An n-point correlation function of the collective monomer density fields 
is given by the expectation value 




(10) 




(11) 



RiiL) = (Ci,(ri)ci2(r2) • • •Q„(r„)) 



1 S^'E 



(12) 



S5/iji(ri) ■ ■ ■5/ij„(r„) 
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or, in more compact notation, 



it:(")(l,...,n) = (c(l)c(2)---cH) 
_ 1 

~E5h{l)---5h{n) ^^^^ 
where c(j) = Ci^(rj). 

Collective cluster functions are related to the correlation functions by a cumu- 
lant expansion. Collective cluster functions are defined, in compact notation, 
as 



5(")(l,...,n)^(c(l)c(2)---c(n))e 

^^Mn^ (14) 
5h{l)---5h{n) ^ ^ 

where the notation (• ■ ■)c denotes a cumulant of the product of variables be- 
tween angle brackets. For example, the two-point cluster function is given 
by the cumulant 'S'(l,2) = (c(l)c(2)) — (c(l))(c(2)), also referred to as the 
structure function. 



2.4 Intramolecular Correlation Functions 
The intramolecular correlation function 

Ma 



^li (r) = { J2 ^o.mh (ri) • • • CamiA^n)) (15) 
m=l 



describes correlations among monomers that are part of the same molecule m 
of species a. In more compact notation. 

Ma 

n(")(l,...,n) = (^c„^(l)---c„^(n)) , (16) 

m=l 



where Cam{j) = Camij(i"j). We also define a function 

Q(")(l,...,n)^^Qi")(l,...,n) . (17) 



The sum in Eq. (17) must be taken only over species of molecule that con- 
tain monomers of all types specified in the argument list i = 
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since Q]^-{r) — otherwise. Note that the single-species correlation function 
Qa(l, . . . , n) and the sum ^('^^(l, . . . , n) defined in Eq. (17) are distinguished 
typographically only by the presence or absence of a species index a. 

The value of . . . ,n) for a specified set of arguments is roughly pro- 

portional to the number density pa = {Ma)/V for molecules of type a, and 
vanishes as pa — 0. To account for this trivial concentration dependence, we 
also define a single-molecule correlation function 

a;W(l,...,n)^QW(l,...,n)/p„ . (18) 



The normalized Fourier transform ci;^"'*(k) may be expressed as an average 
4?(k)= E (e'^^''^-^'""^^'^-)) (19) 



for an arbitrarily chosen molecule m of type a. 

For example, in a homogeneous liquid containing Gaussian homopolymers of 
type a, each containing monomers of type i and statistical segment length 
b, Jajii^) = Kaik'^Nahl/Q), where g{x) = 2[e-^ - 1 + is the Debye 

function. In the limit k — > 0, this quantity approaches ^^^^^^(O) = Nf^. In a 

liquid containing block copolymers of species a, the function tu^ /j(k) describes 
correlations between monomers in blocks i and j. 



2.5 Direct Correlation Function 



We define a collective direct correlation function C(l, 2) for monomers of type i 
and j in a homogeneous fiuid by a generalized Ornstein-Zernicke (OZ) relation 

5-i(l,2) = Q-^(l,2)-C(l,2) , (20) 



or S~j^{k) = VL^^^{k) — Cij(k) in a homogeneous liquid, in which Qr^ is the 
inverse of the two-point intramolecular correlation function QP'^ defined in Eq. 
(17). The random phase approximation (RPA) for 5*"^ in a compressible liquid 
is obtained by approximating —C by the bare potential U . 

The definition of a monomer "type" used here is somewhat fiexible. Because 
all monomers of the same type i arc assumed in Eq. (4) to have the same pair 
potential Uij with all monomers of type j, monomers of the same type must be 
chemically identical. It is sometimes useful, however, to divide sets of chem- 
ically identical monomers into subsets that are treated formally as different 
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types. For example, it may be useful for some purposes to define each "type" 
of monomer to include only monomers of a given chemical type on a specific 
molecular species. If we go further, and take each monomer "type" i to include 
only monomers that occupy a specific site s on a specific species of molecule 
a, then generalized Ornstein-Zernicke equation (20) becomes equivalent to the 
site-site OZ equation used in reference interaction site model (RISM) inte- 
gral equation theories, [10,11,12,13,14,15,16] and C(l,2) to the corresponding 
site-site direct correlation function. 



2.6 Thermodynamic Potentials 



The grand canonical thermodynamic potential $ and Helmholtz free energy 
A are given by 



<l>(^; [h]) = -\nE 
^((M);[/i]) = -lnS + ^/.„(M„) , 



(21) 



Here, square brackets denote functional dependence upon the field h. Corre- 
sponding free energies as functionals of the average monomer densities, rather 
than h, are defined by the Legendre transforms 



r(/.;[(c)])^* + /d(l)Ml)(c(l)) . 
Fi{My,[{c)])^A + Jdil)hil){c{l)) . 

By construction, these quantities have first derivatives (Ma) 
/la — dF/d{M)a, and first functional derivatives 



(22) 
(23) 

—dV/djia and 



5V 



8{c{l)) 



5F 



5{c{l)) 



(24) 



(M) 



Second functional derivatives of T and F are related to the inverse structure 
function by the theorem 



5T 



5{c{l)) 5{c(2)) 



5F 



5{c(l)) 5{c{2)) 



(25) 



(M) 



where S ^(1, 2) is the integral operator inverse of 3^"^^!, 2) 
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3 Field Theoretic Approach 



The field theoretic approach to the statistical mechanics of polymer liquids is 
based upon a mathematical transformation of either the canonical or (here) the 
grand canonical partition function into a functional integral. The functional in- 
tegral representation may be obtained by applying the Stratonovich- Hubbard 
identity [38,39] 

exp |-^c * * c| = J D[J] exp J * C/"^ * J + * c| // (26) 

to the pair interaction energy Uint — * U * c, in which I is the constant 
/ = y^L>[J]exp|-^J*t/-^* j| , (27) 



and 



/ D[J] ^ n / dJiO^) ■ (28) 

ik 



denotes a functional integral of a multicomponent fluctuating chemical poten- 
tial field J= [Ji(r),..., Jc(r)]. 

Inserting Eq. (26) into Eq. (7) for E[h] yields a representation of S as a func- 
tional integral 

S[/i] =/-iyD[J] e^I^''^] (29) 



with 



L[J,h]=lnE[h + zJ]--J*U-^ *J , (30) 

2 



where we have introduced the notation 

E[h] = Tr[ e-^°+~^*'' e"'^^'^ ] (31) 



for the grand partition function of a hypothetical ideal gas of molecules in an 
external field h, with a total energy Uq + Uext[h], but with no pair interaction 
Uint- In Eq- (30), E[h + iJ] is thus the grand-canonical partition function of an 
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ideal gas of molecules in which monomers of type i are subjected to a complex 
chemical potential field hi{r) = /ij(r) + iJj(r), or 

h(l) = h{l)+iJ{l) . (32) 



Hereafter, quantities that are evaluated in such a molecular ideal gas state, 
hke are indicated by a tilde. 

The grand partition function 5[/i] of an ideal gas subjected to a field h is given 
by the exponential 

m = n ( E ^-^^^^^ I = exp (y: x^h]] , (33) 

a [Ma=0 IVl. j y „ y 

where 

Za[h] = J D[R] e-^o[R]+E,^^(R^(^)) (34) 



is the canonical partition function of an isolated molecule of type a in an 

environment in which monomers of type i are subjected to an external field 
hi{r), where Ri(s) is the position of monomer s of type i of the relevant 
molecule. The average molecular number density in such a gas is 

Pa = = — ^ ■ (35) 



Note that InS = J^ai^a), and that the pressure P — InE/V is given by the 
ideal gas law P — Ylia Pa- given by the ideal gas law. 

The value of the intramolecular n-point correlation function Q^^^ in such a 
hypothetical ideal gas will be denoted by Q^"). This quantity may also be 
expressed as a product 

ni'^)(l,...,n)^p„cD(«)(l,...,n) (36) 



where ci;^"''(l, ... ,n) is an ideal-gas single-molecule correlation function. The 
function o)^") is a functional of /i, but is independent of Aq. The function 

0(")(l,...,n) = ^OW(l,...,n) (37) 

a 

is the ideal-gas value of the function fi^^^ defined in Eq. (17). 
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To construct a pcrturbativc expansion of Eq. (29), we will need expressions 
for the functional derivatives of L, and thus of lnS[/i]. The first functional 
derivative of lnS[/i] with respect to h is the average monomer concentration 
field, 



where the expectation value (c(l)) is evaluated in the ideal gas reference state. 
Higher derivatives are given by 



where ^2*^"^ is defined by Eq. (37). 



4 Mean Field Theory 



A simple form of self-consistent field (SCF) theory may be obtained by apply- 
ing a saddle-point approximation to functional integral (29). [19] To identify 
the saddle point of the field J, we require that the functional derivative of Eq. 
(30) for L[J, h] with respect to J vanish. This yields the condition 

iJ^(l) = -|d(2)C/(l,2)(c(2))^ (40) 
or, equivalently, 

h^{l)^h{l)-ld{2)U{l,2){c{2)r . (41) 

Here, and denote saddle point values of the fields J and h, respectively, 
and (c(2))* is the expectation value of c(2) in an ideal gas subjected to a field 
h^. Eq. (41) is the self-consistent field equation for a molecular fiuid in an 
external field h. 

Approximating In S by the saddle point value of L yields a self-consistent field 
approximation for the grand potential $ = — In S as 

^--E{Mar + ll ld{l)d(2)U-\l,2)r(l)r{2) (42) 
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in which {Ma)^ = \aZa[h^]. Combining this with Eq. (41) yields corresponding 
approximations for F and F as 



r = - Y^i^aY + (C) * + U,nt[{c)] (43) 
a 

F = ^(M„) In {^^] + (c) * + C/.„4(c)] , 

where is the saddle-point field corresponding to the specified field (c), and 
U^ntm ^ \ JJd{l)d{2) C/(l,2)(c(l))(c(2)) (44) 

is a mean field approximation for {Umt)- 

In the case of a homogeneous fluid with h = 0, this mean field theory reduces 
to a simple variant of Flory-Huggins theory. In this case, 

M = -UijCj (45) 



where Cj = Y.a ^aiPa is average number density of monomers of type and 
where Uij = jdr' Uij{r' — r). This yields a mean field equation of state 

Pa = exp I - ^ N,Jl,-c, ] , (46) 



and a corresponding Helmholtz free energy density 



^ = i:»^-^) + ',^U,.rc,. (47) 

This is a simple continuum variant of Flory-Huggins theory that (like the 
original lattice theory) neglects all intermolecular correlations in monomer 
density. 



5 Gaussian Field Theory 



Several articles have previously considered a Gaussian approximation to the 
functional integral representation of either the canonical [24,27,28] or grand 
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canonical [35] partition function. A Gaussian approximation for the fluctu- 
ations of J about the saddle-point may be obtained by approximating the 
deviation 



5L[J,^x]=L[J,^x]-L[J\^l] (48) 

of L from its saddle-point value by a harmonic functional 

5L^~j J d{l)d{2) G-\l, 2)5J{1)SJ{2) (49) 

in which 

G-\l,2) = n^^\l,2) + U-\l,2) (50) 



is the second-functional derivative of L about J^, in which D.^'^\l,2) is eval- 
uated at the saddle point. The function G(l, 2) that is obtained by inverting 
is a screened interaction [17,18] analogous to the Debye-Hiickel interaction 
in electrolyte solutions. 



5.1 Grand- Canonical Free Energy 



Evaluating the Gaussian integral with respect to the J field yields an approx- 
imate grand-canonical free energy 

$ = -ln(S) ~ $^ + 5$ (51) 



where is the mean field thermodynamic potential of Eq. (42), and where, 
in a homogeneous fluid. 



d(^> = ^V J lndet[G-^(k)U(k)] 

k 

J lndet[I + 0(k)U(k)] . (52) 

k 

Here, we have introduced the notation 
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for Fourier integrals, and a matrix notation in which bold-face variables indi- 
cate C X C matrices such as G~-^(k), [/jj(k), and the identity — Sij. For 
an inhomogeneous fluid, with h ^ 0, 

S^^-\ndet[G-UU] , (54) 
2 



in which det[...] represents a generalized determinant of a linear integral op- 
erator. We see from the second line of Eq. (52) that 5$ vanishes identically in 
the limit U = of a true ideal gas, in which the saddle-point value already 
yields the correct free energy. 



5.2 HelmhoUz Free Energy 



A Gaussian approximation for the Helmholtz free energy A has been obtained 
in several earlier studies [24,27,28] from the functional integral representation 
of the canonical partition function. Here,we present an alternative derivation 
in which the Gaussian approximation for A is obtained by considering a hy- 
pothetical charging process in which we isothermally turn on the strength of 
the pair interaction while holding the number of molecules fixed. Gonsider a 
model with a re-scaled pair potential energy 

= f / / d{l)d{2) C/(l, 2)c(l)c(2) (55) 



with a charging parameter < a < 1. The Helmholtz free energy A{a) for 
this modified model satisfies an identity 



dA{a) 
da 



1 



(M) 



d{l)d{2) C/(l,2)(c(l)c(2)) 



(56) 



or, in a homogeneous fluid. 



dA{a) 



da 



(M> 



V— __ V f 

—UijCiCj + — J Uij{k)Sij{k;a) 

k 



(57) 



The first term on the RHS of Eq. (57) is the interaction energy of a uniform 
background of monomer density, and the second is a correlation energy. Up to 
this point, the derivation is exact. 

The Gaussian approximation for A in a homogeneous liquid may by obtained 
by approximating Sij{k;a) for for all < o; < 1 by the random phase approx- 
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imation 

S-'(k; a) = [ n'\k) + aU(k) 



(58) 



for the structure function of a system with a rescaled pair potential Q;U(k). 
By substituting approximation (58) into Eq. (57), and integrating the r.h.s. 
of Eq. (57) with respect to a from to 1, we obtain a Helmholtz free energy 
A = A* + SA, in which A* is the Flory-Huggins free energy given in Eq. (47), 
and 

SA = ^V J In det[I + f2(k)U(k)] (59) 

k 

is a Gaussian correction to A. 

Eq. (59) for SA and Eq. (52) for 5$ have the same form, but shghtly different 

meanings: Because 5$(/i) is a grand-canonical free energy, Eq. (52) should be 
evaluated by taking fi*-^-* = Y,aPa^^a \ where pa is a number density obtained 
from the Flory-Huggins saddle-point equation of state for a specified set of 
chemical potentials. Eq. (59) is instead evaluated with Q^^^ = Y.aPa<^^\ in 
which Pa is treated as an input parameter, so that pa may be held constant 
during the "charging" process discussed above. 



6 Perturbation Theory 

The development of a perturbative expansion of $ beyond the Gaussian ap- 
proximation is based upon the expansion of L[J, h] about its value at some 
reference field J* as a sum of the form 

L — Lj-ei + Lharm + La,nh (60) 

where Lj-ef = L[J*,h] is the value of L evaluated at the reference field, and 
-^harm is a harmouic functional of the form 

i^harm ^ -\ j j d{l)d{2) K-\l, 2) 6J{l)6J{2) , (61) 

in which 5J{1) = J(l) — J*(l). The quantity Lanhl"^] is an anharmonic remain- 
der that is defined by Eq. (60). 

The formalism allows us some freedom in the choice of both the reference field 
J* and the kernel K^^. One obvious choice is to take J* to be the saddle- 
point field, and to take = to be the second functional derivative 
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of L[J, h] about its saddle point, as in the Gaussian approximation. Another 
choice will also prove useful, however, and so the values of both J* and K are 
left unspecified for the moment. 

The functional Lanh may be expanded as a functional Taylor series, of the 
form 

^anh = E -T / ^(1) • • • / din) lS")5J(1) • ..SJ{n) , (62) 

n=l n- -f J 



in which 

1 A"7". , r r ^1 

(63) 



J=J* 



Using Eqs. (30) and (60) for L, and Eq. (39) for the derivatives of lnS[/i], one 
finds that 



L«(l) = (c(l)) + jd(2) U-\l,2)tJ*{2) 
L(2)(1, 2) =f^(2)(l, 2) + U-\l, 2) - K{1, 2) (64) 

for n = 1, 2, and 

L(")(l,...,n) (l,...,n) (65) 

for all n > 3. In Eqs. (64) and (65), (c) and fl^^^ are evaluated in the ideal 
gas subjected to a field h — h + iJ*. If J* is taken to be the saddle-point field, 
then = 0. If K'^ is taken to be G'^, then L^^) = 0. 

Given any decomposition of L[J,h] of the form given in Eqs. (60) and (61), 
Eq. (29) for E[h] may be expressed as a product 

'— ' 6 ^ '— 'hajrm'— 'anh ) (^^) 

where 



■^harm 



= C J D[J] e^"^--" (67) 



is a factor arising from Gaussian fluctuations of J in the Gaussian reference 
ensemble. The harmonic contribution corresponds to a free energy 

-lnSha.m=^lndet[X-iC/] , (68) 
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analogous to that obtained in the Gaussian approximation, in which K — G. 
The remaining factor Hanh is given by 



where 



denotes an expectation value in a Gaussian statistical ensemble in which the 
statistical weight for the field J is proportional to e^^'^^^^'^. The inverse of the 
kernel K~^{1, 2) gives the variance 

i^(l,2) = (5J(l)5J(2))ha.m (71) 



of 5 J in this reference ensemble. 

The perturbative expansion of Sanh considered here is based upon a standard 
Wick expansion of Eq. (66) for Sanh as an infinite sum of Fcynman diagrams. 
The required expansion is reviewed briefly in appendix A. The result is an 
expansion Sanh as an infinite sum of integrals. Each integral in this expansion 
may be represented graphically by a diagram of vertices connected by bonds. 
In these diagrams, each vertex that is attached to n bonds represents a factors 
of L^'^\ and each bond represents a factor of —K within the integrand. This 
diagrammatic expansion is discussed in detail in Sec. 8 and several associated 
appendices. 



7 Variants of Perturbation Theory 



In what follows, we consider several variants of perturbation theory that are 
based on the different choices of reference field J* and/or the reference kernel 
K. Throughout this article, we are primarily interested in the calculation of 
properties of a homogeneous liquid, with a vanishing external field, /i = 0. 
Because expressions for correlation functions will be generated by functional 
differentiation of lnS[/i] with respect to /i, however, we must also consider 
states in which h is nonzero, and generally inhomogeneous. Setting the ref- 
erence field J* equal to the saddle-point field J*, for all h, would make J* a 
functional of /i, since the saddle-point fields varies when h is varied. In what 
follows, we instead always take J* to be a spatially homogeneous field that 
is independent of h. The functions L*^") are thus taken to be functionals of a 
field h = h -\- iJ* /m. which J* is independent of h. 
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We consider expansions based upon two different dioices for J*: 

i) Flory-Huggins: We may take J* to be the saddle-point field for a homoge- 
neous liquid with /i = 0, which is given by Eq. (45). This choice, which will 
be referred to as a Flory-Huggins reference field, yields L^^^ — when /i = 0, 
but ^ when h^Q. 

ii) Ideal Gas: Alternatively, we may take J* = 0. This choice, which will be 
referred to as an ideal gas reference field, yields a nonzero value L^^^ = il^^^ 
even when h — 0. 

The functions L and Q that appear in Eqs. (62-63) must be evaluated in an 
ideal gas state with a specified activity for each type of molecule and a 
corresponding molecular number density Pa[h] in this reference state, where 
h — h + J*. The relationship between Aa and Pa thus depends upon our choice 
of reference field J*. In the ideal gas reference state, the molecular number 
density obtained when h = is given by the equation of state pa = \aZa[0]/V 
for an ideal gas in a vanishing potential. In the Flory-Huggins reference state, 
the value of Pa when h — is given by the Flory-Huggins equation of state of 
Eq. (46). 

We also consider expansions based on two different choices for K: 

i) Screened interaction: The obvious choice from the point of view of statistical 
field theory is to take 

K-\l,2)=h{l,2) + U-\l,2) , (72) 

or K = G. This yields a diagrammatic perturbation theory in which the 
propapagator K{1, 2) represents a screened interaction, and in which L^"^^ — 
0. 

ii) Bare Interaction: For some purposes, it is useful to also consider an expan- 
sion in which 



K-\1,2) = U-'{1,2) , (73) 

so that the propagator K = U represents a bare pair potential. This choice 
yields a vanishing harmonic contribution In Sharm = 0, and a nonzero value 

of L(2) =0(2). 



Of the four possible combinations of the above choices for J and two are 
of particular interest: 

Use of a Flory-Huggins reference field and a propagator K — G yields an 
expansion similar to that used in most previous applications of Edwards' ap- 
proach. For a homogenous liquid with h = 0, this choice yields an expansion 
about the Flory-Huggins saddle point with a harmonic contribution In Sharm 
identical that obtained in the Gaussian approximation of Sec. 5. This leads 
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to an expansion of Sanh in which L^^^ and L^^) vanish when h = 0, yielding 
diagrams of —G bonds connecting fl^^^ vertices of all orders n > 3. 

The use of an ideal gas reference field J* — and a propagator K — U yields 
a vanishing harmonic contribution, InSharm = 0, and a perturbative cluster in 
which L(") = for all n > 1. This gives diagrams oi —U bonds connecting 
Qin) vertices of all orders n > 1. This yields a diagrammatic expansion that is 
particularly closely related to the Mayer cluster expansion obtained by Chan- 
dler and coworkers, and that serves as a particularly convenient starting point 
for the derivation of various renormalized expansions. 



8 Diagrammatic Representation 

The integrals that arise in the expansion of S, InS and (in later sections) 
various correlation functions may be conveniently represented by graphs or 
diagrams. Examples of the relevant kinds of diagrams can be found in the 
figures presented throughout the remainder of this article, which the reader is 
encouraged to glance through before reading this section. 

8. 1 Diagrams 

The perturbative cluster diagrams used in this article consist of vertices, free 
root circles, and bonds. A vertex is shown as a large shaded circle with one or 
more smaller circles around its circumference. Each of the small circles around 
the circumference of a vertex is known as a vertex circle. A free root circle is 
a small white circle that is not associated with a vertex. Bonds are lines that 
connect pairs of vertex circles and/or free root circles. 

Each vertex with n associated circles represents a function of n monomer posi- 
tion or wavevector arguments and (generally) n corresponding monomer type 
indices. In this article, vertices usually represent intramolecular correlation 
functions. A vertex with n circles that represents a function i>^"'(l, . . . ,n) is 
referred to as an n-point "t; vertex" . Each of the circles around the perimeter 
of a w vertex may be associated with one coordinate or wavevector argument 
and a corresponding monomer type argument of the associated function. 

Bonds are used to represent functions of two monomer positions, such as a two- 
body interaction ?7y(r — r'), or the Fourier transforms of such functions. Either 
end of every bond must be attached to either a field circle or a free root circle. 
A bond that represents a function &(1,2) of two positions (or wavevectors) 
and two type index arguments is referred to as a "6 bond" . 
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Vertex circles may be either field circles or root circles. Field vertex circles, 
shown as small filled black circles, indicate position or wavevector arguments 
that are integrated over in the corresponding integral. Root vertex circles, 
depicted as small white circles, indicate positions or wavevectors that are 
input parameters, rather than integration variables. In an expansion of an n- 
point correlation or cluster function, each diagram in the expansion must have 
n root circles, each of which is associated with one spatial and one monomer 
type argument of the desired function. 

In perturbative cluster diagrams of the type developed here, each field vertex 
circle must be connected to exactly one bond, and root vertex circles may not 
be connected to bonds. The rules for Mayer cluster diagrams are somewhat 
different, and are discussed in Sec. 22. The use of white and black circles to 
distinguish root and field circles is not actually necessary in perturbative clus- 
ter diagrams, since they can be distinguished by whether they are connected 
to bonds. Some such distinction is necessary in Mayer cluster diagrams, how- 
ever, in which bonds may be connected to both field and root circles. The 
convention is retained here, in part, because it is traditional in Mayer cluster 
diagrams. 

Attachment of one end of a bond to a free root circle is used to indicate 
that the spatial and monomer arguments associated with that bond end are 
a parameters (like the arguments associated with a vertex root circle), rather 
than dummy integration or summation variables. For example, the function 
—U (1, 2) may be represented graphically by a single —U bond attached to two 
free root circles (i.e., two small white circles) with arguments labelled 1 and 
2. Free root circles and vertex root circles are referred to collectively as root 
circles. 

Each diagram F represent a mathematical expression, referred to as the value 
of the diagram, that is given by a ratio 

r = i{T)/s{r) , (74) 

where /(F) is the value of an associated integral, and S{r) is a symmetry 
number. The integral associated with any diagram may be interpreted either 
as a coordinate space integral, which is an integral over the positions associated 
with all of the field vertex sites in the diagram, or as a corresponding Fourier 
integral. In either case, the integrand is obtained by associated a factor of 
v^"\l, . . . ,n) or its Fourier transform with each v vertex with n associated 
field and root vertices, and associating a factor of 6(1, 2) or its transform with 
each b bond. The rules for constructing /(F) for an arbitrary diagram F are 
discussed in appendix B, while the definition and determination of symmetry 
numbers is discussed in appendix C. 
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Physical quantities may often be expressed as sums of the values of all topo- 
logically distinct members of specific infinite sets of diagrams. All references to 
the "sum" of a specified set of diagrams should be understood to refer to the 
sum of the values of all valid, topologically distinct diagrams in the specified 
set. 

The following terminology will be useful to describe sets of diagrams: 

A diagram is connected if every vertex in the diagram is connected to every 
other by at least one continuous path of bonds and vertices. 

A diagram may consist of several disconnected components. A component of 
a diagram a subdiagram whose vertices and free root circles are connected by 
continuous paths only to each other, and not to any other vertices or free root 
circles. A connected diagram thus contains only one component. 



8.2 Grand Partition Function 



A standard analysis of the diagrammatic expansion of functional integral (29), 
which is reviewed in appendices A-C, leads to an expression for Sanh as a sum 



Sanh — 1 + < 



Sum of diagrams of L vertices 
and one or more —K bonds 
with no root circles 



(75) 



The diagrams in this sum need not be connected. In this expansion we asso- 
ciate a factor of —K with each bond, and a factor of L*^") with each vertex 
with n field circles. We include a minus sign in the factor associated with each 
bond as a result of the form chosen for expansion (62), in which a factor of i 
is associated with each Fourier component of 5 J: Each bond thus represents 
the expectation value {lY {5 J {1)5 J {2)) of the product of two factors of i5J. 

The linked cluster theorem [23,40,3] (which is stated in appendix subsection 
F.l) relates In Sanh to the subset of diagrams in the r.h.s. of Eq. (75) that are 
connected: 



InS 



anh 



Sum of connected diagrams of L 
vertices and one or more —K 
bonds with no root circles 



(76) 
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(d) ^^^^-^^ 

OOOO 

(f) (g) (h) (i) 




(j) (k) 



Fig. 1. Examples of diagrams that contribute to expansion (76) of InHanh in di- 
agrams of L vertices and —U bonds. Diagrams (a)-(d) are zero- loop or "tree" 
diagrams, (e)-(i) are one loop diagrams, and (j) and (k) are two loop diagrams. 
Diagrams (f)-(i) are examples of the "ring" diagrams discussed in more detail in 
section 9. Diagrams with the topology of diagrams (b), (d), and the ring diagrams 
(f)-(i) are all excluded from the corresponding expansion of InHanh in diagrams of 
—G bonds and L vertices, because they each contain at least one prohibited L^^) 
field vertex. In an expansion with —G bonds about a Flory-Huggins reference field, 
all diagrams with L^^^ vertices also have vanishing values in the homogeneous state 
with = 0, so that tree diagrams (a-d) and diagram (e) also have vanishing values. 
In this expansion, all tree and one loop diagrams have vanishing values when /i = 0, 
and the only nonvanishing two loop diagrams are those shown above as diagrams 
(j) and (k). 



In a homogeneous fluid, each diagram in this expansion yields a thermody- 
namically extensive contribution to InSanh- 

Expansion (76) is formally applicable to expansions based on any choice of 
reference field J, with any reference kernel K. In expansions with K — G, 
however, in which L^^^ = 0, valid diagrams may not contain field vertices 
with exactly two circles. In expansions about a Flory-Huggins reference field, 
tl^^^ = when h = 0, and so diagrams that contain field vertices with only 
one circle have vanishing values when h = 0. 
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9 Chains and Rings 

Here, we consider the relationship between an expansion of In S about a ref- 
erence field J* with a reference kernel K = U and a corresponding expansion 
about the same reference field with K — G. 

One difference between two such expansions is that InSharm = for any ex- 
pansion with K — U, but In Sharm 7^ when K — U. As a result, expansions 
that use different values of K, which must yield the same value for the total 
InS, must also yield different values for InSanh- 

In addition, the rules for the construction of valid diagrams are different in 
two such expansions: Diagrams containing L field vertices with exactly two 
vertex circles are permitted in a diagrammatic expansion with —U bonds, but 
prohibited in an expansion with —G bonds. In an expansion of InSanh with 
—U bonds, field vertices with two circles can appear only within linear chain 
subdiagrams consisting oi —U bonds alternating with Q^^^ vertices, fike those 
shown in Fig. (2). Furthermore, each such chain subdiagram must either close 
onto itself to form a ring diagram, like those shown in diagrams (f)-(i) of Fig. 
(1), or terminate at both ends at vertices that are not L^^^ field vertices, as in 
diagrams (b) and (d). To show the equivalence of a diagrammatic expansion 
with —U bonds to and a corresponding expansion with —G bonds, we must 
thus consider summations of infinite sets of chain subdiagrams and of ring 
diagrams in the expansion with —U bonds. 

9.1 Chain Subdiagrams 

We first consider a resummation of all possible chain subdiagrams, such as 
those shown in Fig. (2). Let An^ijCk) be the Fourier representation of the sum 
of all chain diagrams that consist of a sequence of n —U bonds connecting n — 1 
^(2) vertices, and that terminate at two free root circles with monomer types 
i and j, with all possible values for the monomer types of the 2(n — 1) field 
circles. The corresponding C x C matrix A„(k) = [74„ jj(k)] may be expressed 
for any n > 1 as a matrix product 



Here, the matrix multipication implicitly accounts for the required summation 
over all possible values of the monomer type indices associated with the field 
circle. 

The screened interaction — Gjj(k) may also be expanded, in the same matrix 
notation, as a geometric series 




(77) 
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Fig. 2. Chain subdiagrams with n = 1,2,3 
- G(k) = -[I + U(k)f2(k)]-^U(k) 

= -U(k) + U(k)^2(k)U(k) • • • 



(78) 



Upon comparing Eqs. (78) and (77), we conclude that G(k) = Z^^i A„(k) or, 
equivalently, that 



G(l,2) 



Sum of chain diagrams of -U bonds 
and O^^'' field vertices with two 
free root circles labelled 1 and 2 



(79) 



By convention, the set of diagrams described in Eq. (79) includes diagram 
Ai{l, 2), which consists of a single —U bond. 

Eq. (79) may be used to expand the integral associated with any valid "base" 
diagram of —G bonds and L field vertices as a sum of integrals associated with 
a corresponding infinite set of "derived" diagrams of —U bonds and L dia- 
grams. Each of the derived diagrams may be obtained from the base diagram 
by replacing each —G bond in the base diagram by any chain subdiagram. To 
turn this geometrical observation into a precise statement about the values of 
sums of infinite set of diagrams, one must must also consider the symmetry 
factors associated with each diagram, and the fact that expansion of all of the 
—G bonds in a base diagram may generate many topologically equivalent di- 
agrams of —U bonds. The required statement (which is proven by considering 
these issues) is given by the bond decoration theorem of appendix subsection 
F.4. In the form required here, it may be stated as follows: 

Theorem: Let F be any valid diagram comprised of any number of L field 
vertices, fl root vertices, —G bonds, and free root circles, which does not 
contain any field vertices with exactly two field circles. Then 



Sum of diagrams that may be derived 
from r by replacing each —G bond in 
F by any chain diagram of —U bonds 
and il^"^^ vertices 



(80) 



To replace a —Gij bond by a chain diagram with free root circles of species i 
and j, we associate the free root circles of the chain diagram with the circles 
that terminate the corresponding —Gij bond in F. 
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To use the above theorem to relate the sum of a set B of base diagrams of 
—G bonds to the sum of a corresponding infinite set A of derived diagrams 
of —U bonds, we must show that every diagram in set A may be obtained by 
expansion of the bonds of a unique base diagram in set B, and, conversely, that 
every diagram that may be obtained by expanding the bonds of a diagram in 
B is in set A. 

It may be shown that every diagram in the set of diagrams that contribute to 
an expansion of In Sanh with K — U may be obtained by derived from a unique 
base diagram in expansion of In Hanh with K — G, except the ring diagrams in 
the —U bond expansion. To see why, note that any diagram of —U bonds with 
no root circles that is not a ring diagram may be reduced to a unique base 
diagram of —G bonds by a graphical procedure in which we erase field vertices 
with two circles, and merge pairs of bonds connected by such vertices, until no 
such vertices remain, and then reinterpret the remaining bonds as —G bonds. 
This procedure fails for ring diagrams, however, because it ultimately yields 
a diagram consisting of a single bond connected at both ends to field circles 
on the same vertex. This diagram cannot be further reduced, and is also not 
a valid base diagram of —G bonds, since it contains a field vertex with two 
circles. The value of In Hanh in an expansion with —G bonds is thus equal to 
the contribution of all diagrams except the ring diagrams to a corresponding 
expansion of In Sanh in terms of diagrams with —U bonds. 

9.2 Ring Diagrams 

Now consider the contribution of the ring diagrams, such as those shown below 
in diagrams (f)-(i) of Fig. (1), to an expansion with K = U . Let i?„ denote 
the sum of all ring diagrams with generic field circles that contain exactly n 
f2(2) vertices and n —U bonds. In a homogeneous fluid, may be expressed 
in a matrix notation as an integral 



in which Tr[- • •] represents the trace of its 2 x 2 matrix argument, and in 
which Sn is a symmetry number. It is shown in appendix section D that the 
appropriate symmetry number for this set of ring diagrams (which is equivalent 
to a single diagram with generic field circles, as defined in the appendix) is 
Sn — 2n for all n > 1. With this symmetry number, we recognize [23] the sum 




(81) 



In Sring =Y,Rn = VY^ 



oo oo 



(-1) 



Tr[(f^(k)U(k))'^] 



(82) 



n=l n=l 



2n 
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as a matrix Taylor expansion of 



In Snng = ^ / Tr ln[I + 0(k)U(k)] 

k 

= ^/lndet[G(k)U(k)] , (83) 

k 

Note that this contribution to In Sanh in the expansion with K = U, for which 
InSharm = 0, is cqual to the harmonic contribution InSharm obtained in the 
expansion with K — G. This, together with the conclusions of the preceding 
subsection, completes the demonstration of the equivalence of expansions of 
In S with K — G and K — U when expanded about the same reference field 
J*. 



10 Collective Cluster Functions 



In this section, we develop expressions for the collective cluster functions by 
calculating functional derivatives of the above expansions of In S with respect 
h. To do so, we must develop graphical rules for evaluating the functional 
derivative of the value of a diagram. We consider separately expansions that 
are obtained by differentiating the expansions of In Sanh ^ ^ sum of diagrams 
oi —U bonds and as a sum of diagrams of —G bonds. 



10.1 Expansion in —U bonds 



The collective cluster function S^^\l, . . . , n) is given in Eq. (12) as a functional 
derivative of In S with respect to h. In an expansion with K — U,\n Sharm = 0, 
so that InS = L[J*, h] + In Sanh [^]- 

Functional differentiation of Eq. (30) for Lref = L[J*,h\ with respect to h at 
fixed J* yields a contribution 



n(")(i,...,n) = 



5"L[J, h] 



5h{l) ■ ■ ■ 5h{n) 



(84) 



j=j* 



for any n > 1. Functional differentiation of Lref thus yields an ideal-gas con- 
tribution to 5(")(l,...,n). 

All further contributions to 5* in this expansion may be obtained by differenti- 
ating values of the diagrams that contribute to expansion (76) of In Sanh The 
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2 2 
1^ . 1^ 



5/z(n+l) "+'1^ 

Fig. 3. Graphical representation of the functional derivative of an L or $7 vertex: 
Differentiation of either an n— point L vertex or and n-point f2 yields an (n+l)-point 
vertex. 

functional derivative with respect to h of the integral /(F) associated with 
an arbitrary diagram F may be expressed as a sum of integral, in which the 
integrand of each integral arises from the differentiation of a factor in the inte- 
grand of /(F) that is associated with single vertex or bond in F. In a diagram 
of L and/or Cl vertices and —U bonds, the factors associated with the vertices 
are functionals of h, while the factors of —U associated with the bonds are 
not. In this case, the functional derivative 6I{T)/6h{l) may thus be expressed 
in this case as a sum of integrals in which the integrand of each is obtained 
by evaluating the functional derivative with respect to h of the factor of L or 
Q associated with one vertex of F. 

The n-point function L*^"-' defined in Eqs. (64) and (65) is equal to (l^"^^ for all 
n > 2 when K = U. The function fi^"^ (1, . . . , n) is an nth functional derivative 
of lnS[/;,] with respect to h = h + iJ . Differentiation of the n-point function 
n^"-* with respect to h at fixed J = J* thus yields a corresponding (n + 1) 
point function 



Moreover, in the cases n = 1 and n = 2 {if K ^ U) in which L^") need not 
equal the difference L*^"^ — fi*^") is always independent of h, and thus does 
not affect the value of the derivative. Consequently, 



for all n > 1. Functional differentiation with respect to h of either an L^"'^ 
vertex or an f^^"-* vertex thus always yields an j^^""*"^) vertex, for all n > 1. 

Differentiation of either an L or f2 vertex may thus be represented graphically 
by the addition of a root circle to the vertex, as shown in Fig. (3), if the result- 
ing vertex is interpreted as an Cl root vertex. Starting from expansion (76) of 
In Sanh; which contains only L field vertices, repeated functional differentiation 
thus generates diagrams in which all of the root vertices (i.e., those to which 
one or more field circles have been added by differentiation) are Cl vertices, 
and all of the remaining field vertices are L vertices. This observation is used 
in appendix subsection F.2 to prove the following theorem 
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Theorem: The functional derivative with respect to h{n) of any diagram F 
comprised of L field vertices, VL root vertices, and —U bonds may be expressed 
as a sum of the values of all topologically distinct diagrams that can be derived 
from r by adding one root circle labelled n to any root or field vertex in F, 
and treating all root vertices in the resulting diagrams as Q vertices and all 
remaining field vertices as L vertices. 

To obtain a diagrammatic expression for S'^^\ we use this theorem to evaluate 
the nth functional derivative of the sum of diagrams in expansion (76) of 
InSanh- By repeatedly differentiating each of the diagrams in this expansion, 
and adding the result to the function 0^"^ obtained by differentiating Lref, we 
obtain an expansion 



n] 



Sum of connected diagrams of 
L field vertices, Vt root vertices, 
and —U bonds, such that each 
diagram contains n root circles 
labelled 1, . . . , n 



(87) 



The root circles of the diagram in the set described above may be distributed 
among any number of root vertices, or may all be on one vertex. The set 
of diagrams described in Eq. (87) includes the diagram consisting of a single 
^(") Yooi vertex and no bonds, which is obtained by differentiation of Lref. The 
statement that the root vertices are labelled 1, . . . , n means that each of the 
root circles in each diagram is associated with a distinct integer j — 1, . . . ,n, 
and that circle j is associated with a position argument or wavevector and 
a corresponding type index argument ij of 5"}"^ (r) or its transform VS^"'^ (k) 



10.2 Expansion in —G Bonds (via Topological Reduction) 

A corresponding expressions for 5*^"-' as a sum of diagrams of L field vertices, 
Q root vertices, and —G bonds may be obtained by resumming chain subdia- 
grams of — C/ bonds and fi^^^ field vertices within the diagrams in expansion 
(87). Every diagram of —U bonds and L vertices that contributes to expansion 
(87) may be related to a unique base diagram of —G bonds. The base diagram 
associated with each diagram of —U bonds may be identified by a procedure 
similar to that discussed in Sec. 9, in which we remove all field vertices with 
two circles, while leaving all other field vertices and all vertices with one or 
more root circles. 

Because the ring diagrams required special care in the expansion of In Sanh, we 
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(a) (b) (c) (d) (e) 



(f) (g) (h) 




(I) 



Fig. 4. Diagrams that contribute to expansion (87) of the two point cluster func- 
tion S'(^)(l,2) in diagrams of —U bonds and root vertices. This expansion is 
obtained by expanding about an ideal gas with J* = using a kernel K = U. 
Diagrams (a) - (e), in which both root circles are on a single vertex, are shown 
in Sec. 11 to contribute to the intramolecular correlation function 0(^''(1,2). The 
infinite sum of diagrams of alternating O^^) field vertices and —U bonds that begins 
with diagrams (a), (f), (g), (h), may be resummed to yield the RPA approximation 
S'-i(l,2) = n-^{l,2) + U{1,2). In the corresponding expansion of 5'(2)(1,2) about 
a Flory-Huggins reference field with a kernel = G, in which bonds represent 
factors of —G, only diagrams (a), (d), (e), (f), (k), and (1) are valid diagrams with 
nonzero values in a homogeneous liquid: In this expansion, diagrams (g) and (h) are 
prohibited because they each contain at least one L^^^ field vertices, while diagrams 
(b), (c), (i), (j) have vanishing values when h = 0, because they each contain at 
least one vertex. 



first consider the one- loop diagrams that are produced by taking nth functional 
derivatives of the sum of ring diagrams in the —U bond expansion of InS. 
Each of the resulting diagrams consists of a ring with a total of n added root 
vertex circles. Each such diagram contain one or more chain subdiagrams of 
alternating —U bonds and Q'^'^^ field vertices. Each such chain subdiagram 
must terminate at both ends at a root vertex with one or more root circles in 
addition to two field circles. Resummation of chain subdiagrams thus yields a 
sum 

f^" In " ■ 

•"'"^ (88) 



6h{l)---6h{n) 
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Sum of diagrams comprised of rings 
of alternating —G bonds and Ct 
root vertices, in which each vertex 
has two field circles and one or more 
root circles, with a total of n root 

circles labelled 1, . . . , n 

Because each remaining vertex contains at least one root circle, in addition to 
two field circles, these (unlike the ring diagrams from which they are derived) 
are valid diagrams of — G bonds. 

Applying the same resummation of chain subdiagrams to all diagrams in the 
—U bond expansion of yields an expression for the collective cluster func- 
tion as 

Sum of connected diagrams of 
L field vertices, Cl root vertices, 
and —G bonds, such that each ) (89) 
diagram contains n root circles 
labelled 1, . . . , n. 



^('^)(l,...,n) = < 



This sum includes the one- loop diagrams described in Eq. (88). Here and here- 
after, it should be understood that a "sum of diagrams is always implicitly 
restricted to valid diagrams, and that valid diagrams of —G bonds cannot 
contain QP'^ field vertices (i.e., field vertices with exactly two vertex circles). 



10.3 Expansion in —G Bonds (via Functional Differentiation) 

Expansion (89) of S'^"'^^ may also be obtained by functional differentiation of 
an expansion of In S as a functional of G*. To calculate derivatives In S in an 
expansion with K — G, we must consider derivatives of Lref, InSharm, and 
InSanh, all of which arc generally nonzero. The functional derivatives of Lref 
are again given by Eq. (84). It is straightforward to show that functional 
differentiation of Eq. (68) for InSharm yields the sum of one-loop diagrams 
described in Eq. (88). 

To evaluate functional derivatives of diagrams — G bonds that contribute to the 
expansion of InSanhj we must take into account that the fact that the factors 
of —G associated with the bonds are functionals of f2*^^)(l, 2; [h]), and are thus 
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Fig. 5. Graphical representation of Eq. (90) for the functional derivative of a —G 
bond. 

also functionals of h. An expression for the functional derivative of —G with 
respect to h may be obtained by differentiating the definition G * G^^ = 6. 
This yields an identity 



6h{3) J J ' ' ' ' ^ ' ' Sh{3) 

d{A) d{5) G{1, 4:)Q^^\3, 4, 5)G(5, 2) 



(90) 



A graphical representation of this identify is shown in Fig. (5). 



In appendix E, we use the graphical rules given in Figs. 3 and 5 to show that 
the functional derivative of any diagram F of L root vertices, Cl field vertices, 
and —G bonds is given by the sum of all diagrams that may be obtained by 
adding one root circle to any vertex or by inserting a Q^^^ vertex with one 
root circle into any bond of F. All of the diagrams that contribute to Eq. 
(89) for S''-"'' may be obtained, using this rule, by functional differentiation of 
Lrcf, InSharm, or InSanh^ Differentiation of Lrcf yields ^2^"^. Differentiation of 
InSharm IS fouud to yield all of the one-loop diagrams described in Eq. (88). 
All of the remaining diagrams in Eq. (89) can be obtained by differentiation 
of one of the diagrams in the —G bond expansion of InSanh- 



11 Intramolecular Correlations 



A diagrammatic expression for the intramolecular cluster function fi^") (1, . . . , n) 
may be obtained by the following thought experiment: We mentally divide the 
Ma molecules of species a within a fluid of interest into two subspecies, such 
that a minority "labelled" subspecies has a number density epa and majority 
"unlabelled" subspecies has a number density (1 — e)pa, with e ^ 1. In addi- 
tion, we imagine dividing the monomers of each chemical type i into labelled 
and unlabelled "types" , by making the monomers on the labelled subspecies of 
molecular species a conceptually distinct from chemically identical monomers 
on any other species of molecule. This mental labelling of particular molecules 
and the monomers they contain is analogous to the procedure used to measure 
intramolecular correlations by neutron scattering, in which a small fraction of 
the molecules of a species are labelled by deuteration. 
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Consider the collective cluster function . . . , n) for a set of labelled 

monomers of types i = {i ...,?„} that exist only on molecules belonging 
to the minority subspecies of species a. In the limit e ^ 0, this function will 
be dominated by contributions that arise from correlations of monomers on 
the same labelled molecule of species a. The single-molecule intramolecular 
cluster function uj^\l, . . . , n) may thus be obtained from the limit 



J:\l,...,n)=]im—S^-\l,...,-e) (91) 



where epa is the number density of the labelled subspecies of a , and S'*^"''(l, . . . ,n 
is the collective cluster function for the specified set of labelled monomer 
types as a function of e. The value of f^a(l, . . . , n) for all of species a, without 
distinguishing labelled and unlabelled subspecies, may then be obtained by 
evaluating the product = Pa<^a- Multiplication of Eq. (91) for cua by pa 
yields 

Qi")(l,...,n) = lini-5(-)(l,...,n;e) (92) 



in which the l.h.s. is the desired value in the unlabelled fluid, with pa molecules 
of type a per unit volume. 

We may use Eq. (87) or (89) to expand the collective cluster function S'*-"^ (1, . . . , 
from the r.h.s. of Eq. (92) as an infinite set of diagrams. In what follows, we 
will refer to monomer types and vertex circles that are associated with the 
labelled subspecies of molecular species a as labelled monomer types and cir- 
cles, and vertices with labelled circles as labelled vertices. Because labelled 
molecules of type a may contain only labelled monomer types, which cannot 
exist on any other species of molecule, the circles associated with each vertex 
must either be all labelled, or all unlabelled. Each labelled Q vertex in a di- 
agram is associated with a factor of epaUJa- Because each labelled vertex thus 
carries a prefactor of e, independent of the number of associated root circles, 
and because the diagram must contain a total of n labelled root circles on one 
or more labelled root vertex, the dominant contribution to S^^\l, . . . , n; e) for 
a set of labelled monomer types in the limit e — > is given by the set of dia- 
grams in which the only labelled vertex in the diagram is a single root vertex 
that contains all n root circles. Dividing the value of any such diagram by e, 
as required by Eq. (92), simply converts the factor of epai^a associated with its 
root vertex into a factor of Qa- The desired expansion of the intramolecular 
correlation function fla for any molecule of species a in the unlabelled fluid is 
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thus 



Sum of connected diagrams of 
L field vertices, —K bonds, and 
a single Cl^ root vertex with 
n root circles labelled 1, . . . , n 



(93) 



Here, the bonds may be either —U or —G bonds, if the corresponding rules for 
the construction of valid diagrams are obeyed. Note that the set described on 
the r.h.s. includes the diagram containing one Cla root vertex and no bonds, 
shown as diagram (a) of Fig. (4). Diagrams (b)-(e) of Fig. (4) also contribute 
to expansion (93), if the root vertex in each diagram is interpreted as an fla 
vertex. 

In what follows, we will also need an expansion of the function Q^'"\l, . . . ,n) 
that is defined in Eq. (17) as a sum of single-species correlation functions. 
Diagrammatic expansions of OS^'^ for different species a but the same set i of 
monomer types contain topologically similar diagrams that differ only as a 
result of the different species labels associated with the Q,a root vertex. The 
sum of a set of such otherwise identical diagrams, summed over all all values 
of species a, is equal to the value of a single diagram in which the root vertex 
is taken to be an Vt vertex, where O*^"^ is defined by Eq. (37). Thus, 



Q(")(l,...,n) = < 



Sum of connected diagrams of L 
field vertices, —K bonds, and a 
single root vertex with n 
root circles labelled 1, . . . , n 



(94) 



The only difference between Eq. (94) and Eq. (93) is the nature of the function 
associated with the root vertex. 



12 Molecular Number Density 



A diagrammatic expansion for the average number density Pa — {Ma) /V for 
molecules of species a may be obtained by applying the identity 



(Ma) = A„ 



(91nS 

dXa 



(95) 
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(b) 



(d) 



Fig. 6. Diagrams that contribute to expansion (97) for pa, using diagrams of —U 

bonds, Q field vertices, and a single Qa root vertex with no root circles. Here, because 
there are no root circles to distinguish root and field vertices, the ^la root vertex is 
distinguished by being whitened. 

to the diagrammatic expansion of InS. Application of the operation Aag|- to 
the function fl vertex yields 



dn(-\l,...,n) ^ d[p,ui''\h ...,n) 

Xn = 2^ Pa 



b 



dpa 



,n) 



(96) 



Here, we have used the fact that the density Pa in the ideal gas reference state, 
which is given by Eq. (35), is directly proportional to at fixed h. The effect 
of the operation upon an 0, vertex is thus to transform it into a Qa root 

vertex. 

By applying this operation to all of the vertices in a diagram, and repeating 
the reasoning used in Sec. 10 to generate equivalent sets of diagrams of —U 
bonds and —G bonds, we obtain a diagrammatic expansion 



Pa^Pa+y{ 



Sum of connected diagrams of 
L field vertices, one or more 
—K bonds, and a single fla 
root vertex, with no root circles 



(97) 



in which (again), —K bonds may be either —U or —G bonds. Examples are 
given in Fig. 6. 



13 Vertex Renormalization 



Wc now discuss a topological reduction that allows infinite sums of diagrams of 
Cl vertices, in which the vertices represent ideal gas intramolecular correlation 
functions, to be reexpressed as corresponding diagrams of fl vertices, in which 
the vertices represent the intramolecular correlation functions obtained in the 
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interacting fluid of interest. The required procedure is closely analogous to 
one described briefly by Chandler and Pratt in their Mayer cluster expansion 
for flexible molecules [8]. It is analogous in a broader sense to the procedure 
used in the Mayer cluster expansion for simple atomic liquids to convert an 
activity expansion into a density expansion [3] . 

The following definitions will be needed in what follows: 

A connecting vertex is one whose removal would divide the component that 
contains that vertex into two or more components. By the removal of a ver- 
tex, we mean a process in which the large circle and associated smaller circles 
representing the vertex and its arguments are all erased. Any bonds that were 
attached to the field circles of the removed vertex must instead be termi- 
nated at free root circles, whose spatial and monomer type arguments become 
parameters in mathematical expression represented by the diagram. 

The components that are created by the removal of a connecting vertex are 
referred to as lohes of the corresponding component of the original diagram. 
A connecting circle thus always connects two or more lobes. A lobe is rooted 

if, before the removal of the connecting vertex, it contained at least one root 
circle (excluding those created by the removal of the vertex, which correspond 
to field circles of the removed vertex), and rootless otherwise. 

An articulation vertex is a connecting vertex that is connected to at least one 
rootless lobe. A nodal vertex is a connecting vertex that connects at least two 
rooted lobes. It is possible for a connecting vertex that is connected to three 
or more lobes to be both a nodal vertex and an articulation vertex. 

In the remainder of this section, we consider a topological reduction of expan- 
sions of various correlation functions and of the molecular density as sums of 
diagrams of —U bonds and f2 vertices. The reduction discussed in this sec- 
tion is not directly applicable to diagrams of —G bonds, for reasons that are 
discussed below. 



13.1 Single- Molecule Properties 

We first consider the reduction of expansion (94) for as an infinite sum 
of diagrams oi —U bonds and VL vertices. Let A be the infinite set of such 
diagrams described in Eq. (94). Let B be the subset of A that contain no 
articulation field vertices, which will be refer to as "base diagrams" . The de- 
sired reduction is based upon the observation that each of the diagrams in set 
A may be derived from a unique base diagram in subset B by "decorating" 
each field vertex with n circles in the base diagram by one of the diagrams 
described in expansion (94) of O^"^ The "decoration" of a vertex is illustrated 
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in Fig. (7): To decorate a vertex v with n circles with a pendant subdiagram 7 
that contains a unique root vertex with n root circles, we superpose the root 
vertex of subdiagram 7 onto vertex v of the base diagram, superpose the root 
circles of 7 onto the circles of v, and blacken any circles that correspond to 
field circles of vertex v in the base diagram. 

This observation about the topology of diagrams is converted into a precise 
statement about their values by the vertex decoration theorem given in ap- 
pendix F.3, which is a generalization of lemma 4 of Hansen and MacDonald. 
As a special case of this theorem, we find the following: 

Theorem: Let F be a diagram of —U bonds and vertices, in which some speci- 
fied set of "target" vertices are Q vertices, none of which may be articulation 
vertices, each of which must contain exactly n field and/or root circles. The 
value of F is equal to the sum of the values of all diagrams that can be ob- 
tained by decorating each of these target vertices with any one of the diagrams 
belonging to the set described on the r.h.s. of Eq. (94) for Q^"). 

To show the equivalence of the sum of a specified set A of diagrams of fl 
vertices with the sum of a corresponding set B of base diagrams, we must show 
that every diagram in A can be obtained by decorating the target vertices of 
a unique base diagram in set B, and. conversely, that every diagram that may 
be obtained by decorating the target vertices of a diagram in B is in A. In 
cases of interest here, the base diagram of a diagram in set A may be identified 
by a graphical process in which we "chp" rootless lobes off of target vertices 
of derived diagrams until the diagram contains no more articulation target 
vertices. 

Applying this reduction to all of the field vertices in all of the diagrams in 



Fig. 7. An example of the "decoration" of a vertex by a pendant subdiagram: Di- 
agram (c) is derived from diagram (a) by decorating the upper field vertex of (a) 
with subdiagram (b). In diagram (c), the n = 5 field vertex to which diagram (b) 
has been attached is an articulation vertex. 




(a) 



(c) 
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expansion (94) for yields the following renormalized expansion 



0(")(l,...,n) = < 



Sum of connected diagrams 
containing any number of Q field 
vertices and —U bonds, and a 
single O root vertex containing n 
root circles labelled 1, . . . , n, with 
no articulation field vertices 



(98) 



This provides a recursive expressions for i^*^"^ vertices for all n. A corresponding 
expansion of Qa may be obtained by replacing the root ^2 vertex by an Cla 
vertex with a specified species a. 



In this reduction, we renormalize the field vertices, creating Q field vertices, but 
leave the root vertex as an unrenormalized Cl vertex. As one way of seeing why 
the root vertex must be left unrenormalized, note that information about the 
intramolecular potential of an isolated molecule enters this set of relationships 
for fi*^") only through our use of an ideal gas correlation function for the root 
vertex: If we had managed to replace the factor of Q associated with the 
root vertex by factor of fl, the resulting theory would not have contained any 
information about the intramolecular potential. 

It may also be shown that the the single root vertex generally cannot be 
taken to be a target vertex of this topological reduction, because a diagram 
of Cl vertices and —K bonds with a single root vertex that is an articulation 
vertex (i.e., is connected to two or more rootless lobes) generally cannot be 
related to a unique base diagram by the clipping process described above. To 
see this, imagine a process in which we first clip all rootless lobes off all of 
articulation field vertices in such a diagram, and then clip off all but one of the 
rootless lobes connected to the root diagram. The resulting diagram contains 
no articulation vertices, but is clearly not a unique base diagram, because the 
process requires an arbitrary choice of which of the lobes connected to the 
root vertex to leave "undipped" in the final step. 

Reasoning essentially identical to that applied above to Q may also be applied 
to expansion (97) of Pa about an ideal gas reference state as a sum of diagrams 
of —U bonds. In this case, taking all of the field vertices to be target vertices 
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yields 



Sum of connected diagrams of fl 
field vertices, one or more —U 



Pa — Pa + y { bonds, and a single root vertex 



of species a with no root circles, 
with no field articulation vertices 



(99) 



13.2 Collective Cluster Functions 



We next consider the topological reduction of expansion (87) for S^'^l In this 
case, we consider the set A of diagrams described in Eq. (87). We divide this 
into a subset in which all n root circles are on a single root vertex, which 
yields expansion (94) for Q, and a set B of all diagrams that contain two or 
more root vertices. Let C be the subset of diagrams in B that contain no 
articulation vertices, which we will refer to a base diagrams. The reader may 
confirm that any diagram in B (i.e., any diagram in A with two or more root 
vertices) may be derived from a unique base diagram in set C by decorating all 
of the vertices of the base diagram, including the root vertices, with one of the 
diagrams in the expansion of fi. The base diagram in set C corresponding to 
any diagram in set B may be unambiguously identified by the clipping process 
described above, which in this case yields a unique result in which none of the 
root or field vertices are connected to rootless lobes. Topologically reducing 
set B thus yields an expression for S'^") as 

,S(")(l,...,n) = Q(")(l,...,n) 

Sum of connected diagrams 

of Q vertices and —U bonds 

with n roots circles labelled (100) 
1, . . . , n on two or more root 
vertices, with no articulation 
vertices 



Here, the contribution f2^")(l, . . . ,n), which may be represented as a diagram 
consisting one vertex and no bonds, is obtained by summing the subset 
of diagrams in A in which all of the root circles are on one vertex. 
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The topological reduction discussed in this section provides a rigorous basis 
for the separation of the calculation of collective correlation functions into: 

(1) A calculation of collective correlation functions, via Eq. (100), as func- 
tionals of the actual intramolecular correlation functions in the fluid of 

interest, and 

(2) A calculation of the effect of non-covalent interactions upon the intramolec- 
ular correlation functions, via Eq. (98). 

In fluids of point particles or rigid molecules, in which the intramolecular corre- 
lations are either trivial or known a priori, the analogous topological reduction 
allows an activity expansion to be converted into a density expansion. In fluids 
of flexible molecules, this conversion is incomplete. In expansion (100) of the 
collective cluster functions, the intramolecular correlation functions Q^") play 
a role partially analogous to that of the number density in the corresponding 
expansion for an atomic liquid. Because expansion (98) of l^i") still contains 
an explicit factor of the molecular activity Aq associated with the root ver- 
tex, however, the reduction does not allow either intramolecular or collective 
correlation functions to be expanded as explicit functions of molecular den- 
sity. Explicit expansions of ^a^\ l^a, and the Helmholtz free energy as 
functions of molecular number density will be derived in a subsequent paper. 

The topological reduction discussed above cannot be directly applied to di- 
agrams of (l vertices and —G bonds, in which the bonds represent screened 
interactions. Such diagrams may contain articulation field vertices for which 
all but two circles are attached by —G bonds to rootless lobes. Clipping all of 
the rootless lobes from such a vertex and re-interpreting the bonds as G bonds 
would yield a field vertex with two circles, creating an invalid base diagram 
of —G bonds. To obtain expansions of the quantities of interest in terms of 
diagrams of Vt vertices with bonds that represent a screened interaction, it 
necessary to first complete the above renormalize of the vertices in an expan- 
sion with —U bonds, and then resum chains of — ^7 bonds to obtain a modified 
screened interaction, in which Q, is replaced by Q. in the definition. This further 
renormalization of the bonds is discussed in Sec. 15. 



14 Leaf Subdiagrams and Shifts of Reference Field 

The renormalization of vertices discussed above was applied to an expansion 
about an ideal gas reference state, with J* = 0, in terms of diagrams with 

—U bonds and VL vertices. In each of the diagrams described in expansion 
(98) of r^i") or in expansion (99) of p^, the factor of Via associated with the 
root vertex must thus be evaluated using the molecular number density Pa[h] 
obtained for a non-interacting gas in a field h = h, rather than the molecular 
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number density obtained from the Flory-Huggins equation of state. 

In botli expansions (98) and (99), tiie root vertex of eacii diagram may be an 
articulation vertex. Among the types of rootless lobes that may be attached to 
the root vertices in these diagrams is a simple "leaf subdiagram, that consists 
of a single —U bond attached at one end to the root vertex and at the other 
to a Q'^'^^ field vertex, as shown in Fig. (8). The terminal Vl^^^ vertex represents 
a factor of monomer concentration, since Q,^j^\r) = (ci(r)). In this section, we 
consider a partial renormalization of the root vertex in each of the diagrams 
of expansion (98) or expansion (99), in which we show that the contributions 
of such leaf subdiagrams may be absorbed into a shift of the value of the field 
h used to evaluate the factor of fl associated with the root vertex. 

The value of a single leaf subdiagram, viewed as a coordinate space diagram 
with generic field circle and one free root circle, is given by the convolution 

Ur)^E [dr'U,,{r-r'){c,{r')) (101) 

The field (p is quite literally a mean field, since it is equal to the ensemble av- 
erage of the potential / d{2) U (1, 2) * c(2) produced by a fiuctuating monomer 
density field. Note that the average monomer concentration (cj(r')) in this 
equation is the actual value in the interacting fluid, because the correspond- 
ing fl^^\2) = (c(2)) vertex has already been renormalized by the reduction 
discussed in the previous section. 

Let r'-"''™^(l, ... ,n) be a diagram with generic fleld circles that consists of a 
single vertex with n root circles, with arguments 1, . . . , n, and m fleld 

circles, to which are attached m pendant leaf subdiagrams, like those shown 
in Fig. (8). The value of such a diagram is given, in compact notation, by an 
integral 

r(n,m) ^ J_ /■ ^(^ + 1) . . . /■ ^(^ + ^) n 02) 

ml J J 

^Q{n+m) _ _ _ ^ ^ + m)(/)(n + 1) • • • 0(n + m) 

where m! is a symmetry number equal to the number of possible permutations 
of the indistinguishable pendant leaves. Let the sum of all such diagrams, with 
any number of leaves, define a function 

oo 

0^")(l,...,n) = ^ r^'^'™) (103) 

m=0 

The function 17'^"+™') (1, . . . , n) in Eq. (102), which is evaluated in an ideal gas 
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1 V'^*¥' 

^ ^ ofB 



Fig. 8. Diagrams consisting of one root fia vertex attached to m leaf sub diagrams, 
which are denoted r("'"^) in the text, for n = 3 and m = 1, . . . , 4. The sum over all 
m, for fixed n, yields the correlation function Oq"^ for a molecule with a specified 
activity Xa in a mean field ^(1) = / d(2)?7(l, 2)(c(2)) 

subjected to a field h — h,is an m-th functional derivative of . . . ,n;h) 

with respect to the applied field h, at fixed molecular activity. The sum ob- 
tained by substituting Eq. (102) into Eq. (103) is a functional Taylor series 
expansion about h — h of a. quantity 

n^:\l,...,n)^n<^:\l,...,n-[h]) (104) 



in which 0*^"^(1, ... ,n;[h]) is an intramolecular correlation function in an ideal 
gas that is subjected to a field 

Ml) = /i(l) + 0(1) . (105) 



The field h is closely related, but not identical, to the saddle-point field used 
in the mean field approximation. The difference arises from the fact that the 
grand-canonical mean field theory uses a self-consistcntly determined estimate 
of (c(l)) for a system with a specified set of chemical potentials, whereas h{l) 
is calculated here using the exact average value (c(l)) to calculate 0. 

In what follows, we will use the corresponding notation 

Pa = p„(A; [h]) (106) 



to denote the value of average molecular number density of species a obtained 
in grand-canonical ensemble for an ideal gas subsected to a mean field h. 
Because of the lack of self-consistency in the definition of h, the values of 
Pa, which are an approximation for the true molecular number densities, are 
generally not consistent with the exact values of (c(l)) used to calculate h. 
In a homogeneous fluid, the only effect of a shift of the reference fleld h used 
to calculate fi^"^ is to change value of the corresponding molecular density 
Pa(A, h) to Pa, without modifying uj^^\ 

It follows from a straightforward application of the vertex decoration theorem 
of appendix subsection F.3, and from our identification of fl as the sum of 
subdiagrams consisting of a root fl vertex with any number of pendant leaves. 
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that 



QH(l,...,n) = < 



Sum of connected diagrams 
containing any number of Q field 
vertices and —U bonds, and a 
single Cla root vertex containing 
n root circles labelled 1, . . . , n, 
with no articulation field 
vertices, and no Q*^^) vertices 



(107) 



Note that the prohibition on il^^^ field vertices in Eq. (107) is equivalent to a 
prohibition on diagrams with leaf subdiagrams. Applying essentially identical 
reasoning to expansion (99) of Pa yields the expansion 



Pa^ Pa + 



V 



Sum of connected diagrams 
containing any number of Q, field 
vertices, one or more —U bonds, 
and a single Cta root vertex with 
no root circles, no articulation field 
vertices, and no OS^^ vertices 



(108) 



In a homogeneous liquid, the 0^ root vertex in Eq. (107) may be replaced by 
a product PocDaj thus absorbing the effect of the mean field (j) into an overall 
pref actor of pa- 



15 Screened Interaction 



In this section, we obtain an expansion in terms of renormalized vertices with 
bonds that represent a redefined screened interaction. Starting from expan- 
sion (98) or (100) of a correlation or cluster function in diagrams with —U 
bonds and VL field vertices, wc may obtain an effective screened interaction by 
summing chains of alternating Q!-"^^ vertices and —U bonds. Repeating the rea- 
soning outlined in Sec. 9, we find that the resummation of all possible chains 
of alternating —U bonds and fi*^^) field vertices yields an effective screened 
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interaction G{1,2), where 



G-\l,2) = n(l,2) + U-\l,2) 



(109) 



The function G is exactly analogous to G, except for the replacement of 
the ideal gas correlation function Jl(^)(l,2) by the corresponding function 
2) = 2) for a molecule in the interacting liquid. 



Resumming chain subdiagrams of —U bonds and Q^^^ field vertices in Eqs. 
(107) and (108) yields the alternate expansions 



Q(-)(l,...,n) 



Sum of connected diagrams of 
—G bonds, fl field vertices, and 
a single Cla root vertex with n 
root circles with arguments 
1, . . . , n, no articulation field 
vertices, and no ^2*^^^ or il^'^^ 
field vertices 



(110) 



and 



Pa 



Pa+y{ 



Sum of connected diagrams of —G 
bonds, fl field vertices, and a single 
fla root vertex with no root circles, 
with no articulation field vertices, 
and no Q^^^ or Q^^^ field vertices 



(111) 
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(a) 





(c) 




Fig. 9. Diagrams of G bonds, field vertices, and a root vertex that that 
contribute to expansion (110) oida\ 

Applying the same reduction to Eq. (100) for the collective cluster functions 
yields 

^(")(l,...,n) = n(")(l,...,n) 

Sum of connected diagrams of Q 
vertices and one or more —G bonds, 

with n root circles labelled 1, . . . , n (112) 

> 

on two or more root vertices, no 
articulation vertices and no Q^^^ 
field vertices 



In Eq. (112), no explicit prohibition on ^2*^^^ vertices is required, because they 
are already excluded by the prohibition on articulation vertices. 

Expansions (110-112) are essentially vertex-renormalized versions of the ex- 
pansions of the corresponding quantities about a Flory-Huggins reference field 
in terms of diagrams of G bonds and Q vertices. Both types of expansion pro- 
hibit diagrams with 1- and 2-point field vertices, and thus yield diagrams with 
similar topology, in which the vertices and bonds have similar meanings. They 
differ only as a result of the replacement in the renormalized expansions of VL 
field vertices by Vt vertices, and of G bonds by G bonds, and the corresponding 
prohibition in the renormalized expansions on diagrams with articulation field 
vertices. 



16 Bond-Irreducible Diagrams 



The diagrammatic expansion of the two-point cluster function 5'(1, 2) may be 
simplified by relating 5 to a quantity A(l, 2) that wc define as the sum of 
all bond-irreducible two-point diagrams. A diagram is said here to be hond- 
irreducible if it is connected and cannot be divided into two or more compo- 
nents by removal of any single bond. By the removal of a bond, we mean the 
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erasure of the bond, the erasure of any free root circles to which it is attached, 
and the transformation of any vertex field circles to which it is attached into 
vertex root circles. 



We define 



Sum of bond irreducible diagrams 
of 1 or 2 Q root vertices and any 
A(l, 2) = < number of L field vertices and —K 
bonds, with 2 root circles labelled 
1 and 2 



(113) 



The set of diagrams described in on the r.h.s. of Eq. (113) also includes a 
subset of diagrams in which both root circles are on a single vertex. Among 
these is the trivial diagram consisting of a single 0*^^^ vertex and no bonds. 
The sum of this infinite subset of diagrams yields the diagrammatic expansion 
of Q^^^(l, 2). We may thus distinguish inter- and intra-molecular contributions 
to A(l, 2) by writing 

A(l,2) =0(2)(1^2) + E(1,2) (114) 



where 



Sum of bond irreducible diagrams 
with two VL root vertices, any number 
2) = of L field vertices and —K bonds, 
with 2 vertex root circles labelled 1 
and 2 on different root vertices 



(115) 



Renormalizing both vertices and bonds yields the equivalent expansion 



E(l,2) 



Sum of bond-irreducible diagrams 
of Vt vertices and —G bonds, with 
2 root circles labelled 1 and 2 on 
different vertices, with no 
articulation vertices 



(116) 
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(a) (b) (c) 



Fig. 10. Examples of diagrams of G bonds and f2 field vertices that contribute to 
expansion (116) of S(l,2). 

in which field vertices with two circles are prohibited. 

The expansion of S{1, 2) in terms of Vt vertices and —U bonds may be ex- 
pressed diagrammatically as an infinite set of chains of alternating —U bonds 
and A field vertices, terminated at both ends by A root vertices. This may be 
expressed algebraically as an infinite series 

5 = A-A*C/*A + --- 



Summing the series, or solving the recursion relation, yields an expression for 
the inverse structure function as 



Eq. (118) is an generalization of the random phase approximation for in 
which VL is replaced by A. The generalized Ornstcin-Zcrnick relation of Eq. 
(20) is a different generalization in which U is instead replaced by — C. 



17 Direct Correlation Function 

We now construct a diagrammatic expansion of the Ornstein-Zernicke direct 
correlation function C defined in Eq. (20). This quantity may be related to 
the function A by equating the r.h.s.'s of Eqs. (118) and (20) for S'"^(l,2). 
This yields 



= A-A*C/*5 . 



(117) 



5-^(1,2) = A-^ 



(1,2) + C/(1,2) 



(118) 



C(l,2) 



t/(l,2) + AC(l,2) 



(119) 



where 



AC(l,2) = -A-^(l,2)+l^-i(l,2) , 



(120) 



and where A(l, 2) = ^(1, 2) + E(l, 2). 
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By expanding A ^ = [Q + E] -"^ as a geometrical series, AC(1, 2) may also be 
expressed as a convolution 



(121) 



in which 



r=E-E*Q~^*EH 

= E-E*Q-^*T . 

Summing this series, or solving the recursion relation, yields 

r-^(l,2) = E-^(l,2) + Q-i(l,2) 



(122) 



(123) 



Eq. (122) provides a useful starting point for the construction of an explicit 
diagrammatic expansion of T. 



17.1 Diagrammatic Expansions 



To express Eq. (122) as a sum of diagrams, we depict each factor of —VL~^ in 
the infinite series as a — bond, which will be represented by a dotted line. 
This yields an expansion of T as 



Sum of diagrams containing any 
number of Q vertices, —G bonds 
bonds, and — Q"-*^ bonds, with two 
root circles with arguments 1 and 
and 2 on different vertices, with 
no articulation vertices, such that 
the diagram cannot be divided into 
components by cutting any one —G 
bond, but can be so divided by 
cutting any one —VL~^ bond 



r(i,2) 



(124) 



As with any expansion in terms of —G bonds, field vertices with exactly two 
circles are implicitly prohibited. The above description applies equally well 
to an expansion of T(l,2) in terms of diagrams of —U bonds, rather than 
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(a) (b) 
(c) 

o--oo--<I> 

(d) 

Fig. 11. Examples of diagrams that contribute to the expansion (124) for T(l,2). 
Here, vertices are f2 vertices, sohd hnes are —G bonds, and dashed hnes are —Vt~^ 
bonds. A corresponding expansion of AC(1, 2) may be obtained by attaching exter- 
nal — r2~^ bonds to both root sites of each diagram in this expansion. 

— G bonds, if Q*^^^ field vertices are allowed. Corresponding expansions of Eq. 
(121) for AC(1, 2) may then be obtained by adding two external —Vt~^ bonds 
to each diagram in the expansion of T, as shown in Fig. (11). 

An alternative expansion of T may be obtained by viewing each diagram in 
expansion (124) as a necklace of "pearls", where a pearl is a bond-irreducible 
subdiagram of —G or —U bonds and Vl vertices that contains no nodal vertices. 
The pearls of a diagram in expansion (124) are thus the disjoint components 
that would created by removing all — bonds and all nodal vertices. For 
example, in Fig. (11), diagrams (b) and (c) each contain two identical pearls, 
which are connected by a nodal vertex in diagram (b) and by an — Q"^ bond 
in diagram (c). Diagram (d) contains 4 pearls. 

Each diagram in Eq. (124) for T(l,2) may thus be viewed as a string of 
pearls, in which each pair of neighboring pearls may be connected by either 
a nodal vertex or by an — bond connecting two Vt vertices. The set of all 
such diagrams can be generated from the set of bond-irreducible two-point 
diagrams that contributes to Eq. (116) for E by replacing every nodal vertex 
in any diagram in expansion (116) by either the original Vl vertex or by a 
subdiagram consisting of two vertices connected by a — Q bond. Let 

N(n,n') = Q(n+n') _^ j-[(n,n') ^^25) 

denote the sum of the two such subdiagrams shown, for n = n' = 3, in 
Fig. (12). Here, n*^"'"') denotes the value of a subdiagram consisting of an 

vertex and an Q^" vertex connected by a — bond, in which one 
vertex contains n root circles labelled 1, . . . , n and the contains 
n' root circles labelled 1', . . . , n', as shown for n = n' = 3 in the second 
diagram of Fig. (12). Let a nodal N vertex be one in which the function 
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, n') is associated with an (n+n')-point nodal vertex that 
is connected to one pearl by bonds connected to n field circles 1, . . . , n and to 
a second through n' field circles labelled 1', . . . , n' . An alternate expansion of 
T may be obtained by replacing all of the nodal Q vertices in expansion (116) 
of E by nodal N vertices. The corresponding expansion of AC is obtained 
by adding two external — fi"^ bonds to each diagram in the expansion. For 
example, replacing the nodal Vt vertex in diagram (b) of (10) with an nodal 
vertex and then adding two external — fi"-*^ bonds yields a diagram whose 
value is equal to the sum of the values of diagrams (b) and (c) in Fig. (11). 
Geometrical series (122) for T(l, 2) may thus be expressed as a sum 

Sum of bond-irreducible diagrams of 
vertices and —G bonds, with two root 
circles labelled 1 and 2 on different 
T(l,2) = < n vertices, in which all nodal vertices ) (126) 
are N vertices and all other vertices 
are Q vertices, with no articulation 
vertices 



The above description applies equally well to a corresponding expansion of 
T in terms of diagrams oi —U bonds, aside from the removal of the implicit 
prohibition on Vt^"^^ field vertices. 



17.2 Atomic Limit 



Expansion (126) may be used to recover a simpler expression for C in the limit 
of a hquid of point particles. For this purpose, it is convenient to consider the 
expansion of T in terms of —U bonds, rather than —G bonds. In a homoge- 
neous mixture of point particles, the normalized Fourier tranform fi-"^(k) is 
nonzero only when ii = 12 = . . . = in = i and ki + ■ ■ ■ + k„ = 0, and in this 
case is equal to the density pi of particles of type i. Using these relations, it is 
straightforward to show that the function A'^"'"') defined in Eq. (125) vanishes 
identically in the atomic limit. 

Because the function A^ that we associate with nodal vertices in (126) for T 
vanishes, this expansion (126) for T thus reduces in the atomic limit to a sum 
of bond-irreducible diagrams of —p vertices and either —G or —U bonds with 
two root circles on different root circles, with no articulation vertices or nodal 
vertices. Such diagrams can thus have no connecting vertices. Adding two 
external legs with values J7^-^(k) = Sijp^^ to each diagram in the expansion of 
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2' + 



I*" 



Fig. 12. Graphical representation of the quantity N defined in Eq. (125). The quan- 
tity shown is N^^'^\l, 2, 3, 1', 2', 3'), where vertices are Q^") vertices, and the dashed 
line represents a —Q"^ bond. The second diagram is E^^'^^ 

T yields a corresponding expansion for AC. The factors of p^^ associated with 
these external legs cancel the factors of pi associated with the root vertices, 
thus converting the root vertices into "1 vertices" that introduce a factor of 
unity into the integrand of the corresponding integral. Such root 1 vertices 
are represented as white circles in the conventional representation of cluster 
diagrams for atomic liquids. [3] By considering the resulting expansion of AC 
in the atomic limit in terms of diagrams of —U bonds, and adding to this the 
diagram consisting of a single —U bond to obtain C = —U + AC, we obtain 
an expansion 



C(l,2) 



Sum of diagrams of one or more —U 
bonds, any number of field p vertices, 
and two root 1 vertices labelled 1 
and 2, with no connecting vertices 



(127) 



for liquids of point particles. Note that the set described on the r.h.s. includes 
the diagram consisting of a single —U bond, since this description does require 
that the diagrams be bond-irreducible, but the prohibition on connecting ver- 
tices automatically excludes all bond-reducible diagrams other than this trivial 
one. Eq. (127) is a perturbative cluster diagram variant of a well known ex- 
pansion for C(l,2) as a sum of Mayer cluster diagrams, given in Eq. (5.2.16) 
of Hansen and MacDonald [3]), in which the Mayer diagrams are required 
to satisfy the same topological constraints as those specified for perturbative 
cluster diagrams in Eq. (127). 



18 Long Range Interactions and the Loop Expansion 



The field-theoretic approach followed here leads naturally to a loop expansion 
about a mean field saddle point. The resulting expansion is useful, however, 
only if fluctuations about the mean field solution arc small. In fluids of point 
particles, the mean field theory of Sec. 4 is known [41,42,43] to be exact in the 
so-called Kac limit of long-range interactions and/or high particle densities, 
in which each particle interacts weakly with many other particles. 
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To define this limit, we consider a class of models in which the pair interaction 
between monomers or (in an atomic liquid) atoms is of the form 

U,,{v-v')=U,,^'F{^\v-v'\) , (128) 



Here, Uij is an interaction strength with units of volume, 7"^ is an adjustable 
range of interaction, and -F(x) = -F(|x|) is a dimensionless function that is 
required to satisfy a normalization condition / dx F(x) = 1. This parameter- 
ization yields a Fourier transform 

C/,,(k) = C7,,-F(k/7) , (129) 



in which F(k) = / dx e*^'^F(x) is the Fourier transform of -^(x), and in 
which F(0) = 1 as a result of the normalization assumed for F(x). To further 
simplify the argument, we will assume in what follows that F(x) is a function 
that falls off smoothly and rapidly with dimensionless separation |x|, such 
as a Gaussian, and that, correspondingly, F(q) decreases rapidly to zero for 
|q| > 1. 

The mean field theory for simple atomic fluids with interactions is known 
[41,42,43] to be asymptotically exact in the limit 7 ^ of an infinitely long 
range, infinitely weak potential. In this limit, the first correction to mean field 
theory in an expansion in powers of 7 is given by the Gaussian approximation 
discussed in Sec. 5. We show below that this 7-expansion corresponds to a 
loop expansion of the perturbative cluster expansion. Interestingly, the loop 
expansion may be shown to remain valid even for nearly incompressible fiuids 
with very strong repulsive interactions, if the pair potential is of sufficiently 
long range. 



18.1 Nearly Incompressible Atomic Mixtures 



As a simple example that retains some of the physics relevant to a dense 
molecular liquid, we first consider a simple model of a nearly incompressible 
binary mixture of point particles. We consider a mixture with a total particle 
number density c = ci + C2, with interactions of the form given in Eqs. (128) 
and (129). The interaction matrix Uij is taken to be of the form 



(130) 
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This yields an interaction that reduces in the hmit of slow spatial variations 
g -C 7 to a continuum approximation 

Uint jdv |^5o(ci + c^f + X0C1C2I . (131) 



In this form, it is clear that is proportional to a mean-field compression 
modulus and that Xo is a measure of the mean field enthalpy of mixing. The 
quantities Bq and xo have units of volume, and may be used to define di- 
mensionless quantities Bq = Bqc and Xo = Xoc The quantity Xo is a "bare" 
dimensionless Flory-Huggins parameter. 

Consider the screened interaction G for this model in the limit i?o ^ 1 of a 
nearly incompressible liquid. In the limit /c <^ 7, for which F(k/7) ~ 1, taking 
the limit Bq ^ 00 of an infinitely strong bare repulsion yields a screened 
interaction 



G(k) 



c - 2xoCiC2 



1 - 2xoC2 1 
1 1 - 2XoCi 



(132) 



For any large but finite Bq, however, we will assume that a rapid decrease in 
F(k/7) for k ^ 7 causes both the bare interaction U and the corresponding 
screened interaction G to become negligible for |k| 3> 7. In what follows, 
to characterize the relative magnitude of different contributions to the free 
energy at the level of power counting, this behavior will approximated by a 
sharp cutoff in which G(k) is given by Eq. (132) for all < 7, and G(k) = 
for A; > 7. In an interaction site model, the inverse range of interaction 7 thus 
acts as a cutoff wavenumber for the screened interaction. 

The Gaussian contribution to the Helmholtz free energy in the limit i? 3> 1 is 
given in this approximation by an integral 

j ln[(c - 2XoCiC2)5o] . (133) 

|k|<7 



This yields a contribution to the total free energy per particle of order (7^/c) ln(S), 
and a contribution to the free energy of mixing of order Xo7'^ per particle. For 
stable liquids, with Xo ^ 1- the Gaussian contribution to the free energy of 
mixing is thus small compared to the mean field contribution whenever 7^ -C c. 
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18.2 Loop Expansion 



We now consider the loop expansion of In S and its functional derivatives for 
this model. Consider the contribution of an arbitrary connected diagram F of 
—G bonds and c vertices to an expansion of In S about the mean field reference 
state. We consider the Fourier space interpretation of T/V, which represents 
a contribution to the pressure P = InS/V. Let the diagram of interest con- 
tain hb G bonds and ny vertices. In a homogeneous fluid of point particles, 
each vertex simply represents a factor of c, independent of the number of 
field and/or root circles attached to the vertex. The number of indepen- 
dent wavevector integrals in the Fourier representation of the diagram in a 
homogeneous liquid, or the number of loops, is [23,44] 

riL = riB - ny + 1 (134) 



To characterize the dependence of the value of the diagram upon 7 and c, at the 
level of power counting, we consider an approximation with sharp wavevector 
cutoff, in which we use Eq. (132) for — G(k) for bonds with |k| < 7, and set G — 
for all |k| > 7. In this approximation, in a system of point particles, the only 
length scale in the integrand is cutoff length 7"^, and so it is convenient to non- 
dimensionalize all wavevectors by 7. After non-dimensionalization, the value 
r/V will contain a pref actor of order c from each vertex, a pref actor of 1/c from 
each —G bond, and a prefactor of 7^ from the measure of each independent 
wavevector integral (or loop), which multiply a non-dimensionalized Fourier 
integral. Using Eq. (134), we find an overall prefactor of order 

I'b'/cr-' (135) 



for any diagram with loops, where c is an average particle concentration. 
Here, for the purposes of power counting, we have neglected the distinction 
between c and the Flory-Huggins result for c. This prefactor multiplies a di- 
mensionless Fourier integral in which the dimcnsionlcss wavevector k/7 asso- 
ciated with each bond of wavevector k is restricted to values |k|/7 < 1, and in 
which the integrand is a product of factors of Gjj(k)/c that can be expressed 
as a function of Xo, (pi, and 02, where 0j = q/c. The Gaussian approxima- 
tion to the free energy density, which involves a single wavevector integral, 
has a prefactor of 7^, and thus conforms to the rule for a one-loop diagram. 
Eq. (135) implies that the loop expansion will rapidly converge if and only if 
C7~^ 3> 1. This is equivalent to the requirement that many particles lie within 
a distance 7"^ of any test particle. 
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19 Coarse-Grained Models of Dense Polymer Mixtures 



In previous studies, the field theoretic approach has often been used to study 
coarse-grained models of dense multi-component polymer liquids. To clarify 
the assumptions underlying this application, we now consider a class of mod- 
els in which each "monomer" represents a subchain of g chemical monomers, 
and in which such coarse-grained monomers interact via an effective two body 
interaction. For g ^ 1, these monomers are diffuse, strongly overlapping ob- 
jects, and presumably exhibit an effective interaction with a range of inter- 
action comparable to the coil size y^a of a corresponding subchain, where a 
is the statistical segment length of a chemical monomer. The number density 
of coarse-grained monomers in such a model c = l/{vg), where u is a vol- 
ume per chemical monomer. The number of coarse-grained monomers that lie 
within a range of interaction 7"^ ~ ^/ga of any one monomers is comparable 
to the number of subchains that interpenetrate any subchain of length g. This 
quantity is proportional to the ratio 07"^ = h/p, in which 

p = v/a^ (136) 



is the packing length of the melt and h = y/ga is the coarse-grained statistical 
segment length. 

For very high molecular weight polymers, one can envision a model with 
coarse-grained monomers of size 1 <^ g <^ N in which the monomers are 
large enough to overlap, but small enough to allow the description of compo- 
sition fluctuation with wavelengths of order the coil size y/N a or longer. Such 
a model is a potentially useful starting point for studying universal aspects of 
polymer thermodynamics that lie beyond mean field theory, such as the effects 
of fluctuations near the critical point of a polymer blend or the order- disorder 
transition of a block copolymer melt. 

A mean field approximation for the total free energy of such a coarse-grained 
model might be expected to be valid in the limit 7^^ ~ 6 ^ p in which the 
coarse-grained monomers strongly overlap. This assertion is based upon the 
idea that the total free energy of a liquid of polymers consisting of monomers 
depends upon the range of interactions in a manner qualitatively similar to 
that of a corresponding liquid of point particles, and that this criterion corre- 
sponds to the citerion cj~^ ^ 1 found for an atomic liquid. 

To justify this statement more carefully, we now consider a loop expansion 

of InH for consider a model containing coarse-grained monomers connected 
by Gaussian springs of statistical segment length b, with a coarse-grained 
monomer number density c. We asssume a pair interaction of the form given 
in Eq. (130), and consider the nearly incompressible limit Bqc ^ 1. Rather 
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than immediately requiring that 7"^ ~ b, as suggested by the above physical 
arguments, we consider separately the limiting cases 7^^ ^ b and b ^ 7^^ 
separately, and show that both limits extrapolate to the same results for a 
model with ~ 6 ^ p. 

In the limit b ^ 7^^, polymers may be adequately described as continuous 
Gaussian threads over lengths of order the range of interaction 7^^. In this 
limit, we may use the Gaussian thread model for intramolecular correlations 
functions in order to analyze a loop expansion of In S at the level of power 
counting. We consider the limit of constant monomer density c and infinite 
degree of polymerization. In this limit, it may be shown for a homopoly- 
mer a with monomers of statistical segment length b that i7^"-'(k)/c is a di- 
mensionless function of {ki6, . . . , k„6}, for which f2*^")(Aki, Ak2, . . . , Ak„) = 
A2("-^)Q(")(ki, . . . , k„)/c. The magnitude Q(")(k) is thus of order 



when all of its wavevevector arguments have magnitudes of order k. For the 
familiar case n — 2, the Debye function yields Q^^\k) ~ 12c/(/c^6^) in this 
high wavenumber limit. In the same limit, the screened interaction (j(k) may 
be shown [35] to be of order 

G{k) ~ {kbf/c (138) 



for all -C /c <^ 7 for systems with Bq ^ 1 and Xo-^ ^ 1) where R ~ \/Nb. 
Consider a diagram with ub —G bonds and ny ^ vertices in the expansion 
of InSanh- For this purpose, it is convenient to associate each of the n factors 
of {kby in the denominator of the r.h.s. of Eq. (137) for the factor of Jl*^"^ 
associated with a vertex with one of the n attached bonds. This leaves an 
overall factor of order ck^b^ for each vertex, independent of n, and a factor 
of order l/{ck%'^) for each bond, which is obtained by dividing Eq. (138) for 
G by two factors of {khY that are borrowed from the two attached vertices. 
Using these estimates, it may be shown by a straightforward power-counting 
argument similar to that given above for an atomic liquid that if the integral 
associated with a generic connected ni-loo]) diagram with no root circles con- 
verges in the limit A?" ^ 0, it is dominated by wavevectors of order 7, and that 
the integral is of order 

7^(7P)"^"^ ■ (139) 



In this string-like limit, we thus expect the loop expansion of InS to rapidly 
converge whenever 7p -C 1, i.e., whenever the packing length is much less 
than the range of interaction, so that each point on a chain interacts directly 
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with many other chains. In the physically relevant limit 7"^ ~ b, in which 
p = l/{dP') ~ Eq. (139) reduces to Eq. (135) for the magnitude of 

nx,-loop contributions in a fluid of disconnected monomers. 

We next consider the opposite limit, 6 ^ 7"^, in which consecutive monomers 
within a chain are connected only very loosely, by harmonic tethers with a 
characteristic bond length h much longer than either the typical distance c^^^^ 
between neighboring monomers or their range of interaction 7"^. In this limit, 
the connectivity of the polymers introduces only a small perturbation of the 
structure of a reference liquid of disconnected monomers with the same den- 
sity and the same pair interaction. The introduction of such loose bonds is ex- 
pected to signficantly affect only long-wavelength correlations, at wavelengths 
of order h and greater. Because corrections to the mean field free energy are 
controlled primarily by correlations over distances of order 7"^ and less, the 
free energy of such a polymer liquid should be almost identical to that of 
corresponding monomer liquid, and should thus be describable be the loop 
expansion constructed above for an atomic mixture, in which the magnitude 
of a generic riL-loop contribution to InS is given by Eq. (135), in which c is 
the monomer density. This expression may thus be be obtained for a system 
with h ~ 7"^ ^ p by extrapolating from either the limit h <S 7"^ or from 

6> 7~^- 



20 One Loop Approximation for Two-Point Correlations 

As an example of the formalism, and of its relationship to earlier field-theoretic 
studies, we now derive expansions of r2jj(k), a;jj(k), and Cjj(k) for a nearly 
incompressible polymer liquid to first order in a renormalized loop expansion, 
using one-loop diagrams of Q vertices and G bonds. The calculation given here 
is limited to a class of systems with two types of monomer, which includes 
both binary homopolymer blends and diblock copolymer melts. 



20.1 Incompressible Limit 

We are interested here primarily in effects of composition fluctuations in nearly 
incompressible liquids. We thus consider a model with a monomer-monomer 
pair potential of the form given in Eqs. (128-131), and assume that Bqc » 1, 
where c — ci-\-C2 is the total monomer number density. 

By following the reasoning applied above to a nearly incompressible atomic 
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mixture, we obtain a screened interaction 



G(k) ~ 



0+(k)-2xo|^^(k)| 



1- 2x0^22 (k) 1 + 2x0^^12 (k) 
1 + 2x0^^21 (k) 1- 2x0^^11 (k) 



(140) 



for Bqc » 1 and /c <^ 7 ^, and G'(k) ~ for A; » 7 and large but finite values 
of B. Here, 



Q+(k) = Qn(k) + Q22(k) + Qi2(k) + Q2i(k) 
|Q(k)| = Qn(k)Q22(k) -Qi2(k)Q2i(k) 



(141) 



are the sum of elements of f^ij(k) and its determinant, respectively. 

Ornstein-Zernicke equation (20) gives the structure function >S'jj(k) as the 
inverse of a matrix n^j^(k) + Uij + ACij(k). In the incompressible limit Bqc — > 
00, this matrix inverse approaches 



S(k) = 



\m\ 



l]+(k)-2xa(k)|l](k)| 



+1 -1 

-1 +1 



(142) 



where 

Xa(k) = xo - ACi2(k) + ^[ACn(k) + ACn(k)] (143) 

is the wavenumber-dependent apparent x parameter identified previously by 
Schweizer and Curro. [13] 



20.2 Intramolecular Correlations 



The only diagrams that contribute to Eq. (110) for Q;j^(k) to first order in 
a loop expansion are the trival tree diagram consisting of a single vertex, 

shown as diagram (a) of Fig. 9, and the one-loop diagram shown as diagram 
(b) in the same figure. In the renormalized expansion considered here, the one 
loop diagram is evaluated using an Cla = Pa^a vertex and a single —G bond. 
The value of this one- loop diagram, in a Fourier representation, is Vpalif^\ 
where 

/Sf (k) -k, q, -q)G«(q) ■ (144) 

q 
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Summation over all monomer types k and I that exist on molecules of type a 
is implicit. The sum of the tree and one- loop diagram yields a 1-loop approx- 
imation 

Q2^(k)^p„[a;15,.(k) + 7£f(k)] (145) 



where pa is a molecular density that is related to Xa by the Flory-Huggins 
equation of state. 

To calculate the corresponding single- molecule correlation function LVa,ijO<-), 
we divide the above approximation for Qa,ijO^) by a corresponding one- loop 
approximation for the molecular density pa at a specified chemical potential. 
The only 1-loop contribution to Eq. (Ill) for the ratio pa is one with the 
topology of that shown as diagram (c) of Fig. (6). In the renormalized expan- 
sion used here, the bond in this diagram is taken to be a G bond, and the root 
vertex to be a = PaOJa vertex. This yields a 1-loop expansion 

p„~p4l + /f ] (146) 



in which 



4°^ = -^/'^S(q,-q)G«(q) (147) 



To obtain an approximation for uj^Jj, we take the ratio of Eqs. (145) (146), to 
obtain 

- ^"^W + . (148) 

1- -\- la 



To obtain a one- loop approximation, we then expand [1 -|- as a geomet- 
rical series in J^°\ and truncate after the first term in the series. This yields 
an expression 



in which the last two terms each contain an integral with respect to a single 
wavevector, thus yielding a one-loop approximation. The last two terms in Eq. 
(149) may be combined into a single integral 

45,.(k) ulS^ik) -\j V;l^.,,(k, -k, q, -q)G,,(q) , (150) 

q 
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in which 



i^ajjkii^, -k, q, -q) = oj^^IjkiiK -k, q, -q) - oj^SjO^^ -^)^aU<i, -q) ■ 

It is straightforward to show, using Eq. (18) for ciJ^"^(k), that i^tljkii^i q> ~q) — 
0, and thus that this approximation satisfies the exact result 




Eq. (150) has been obtained previously by Barrat and Predrickson [26] for the 
case of a diblock copolymer melt. 

20.3 Direct Correlation Function 

The direct correlation function is related by Eq. (121) to the quantity Tjj(k). 
Expansion (126) of Tjj(k) contains diagrams that have the same topology as 
those in expansion (116) for Eij(k), but that generally have different values as 
a result of the use of N nodal vertices in Eq. (126), rather than Vt vertices. The 
only one-loop diagram in the expansion of cither T or S is diagram (c) Fig. 
(10). This diagram contains no nodal vertices, and thus has the same value in 
either expansion. This diagram yields a one-loop approximation for either T 
or E as an integral 

T,,.(k) ^ i I l]S2(k, q, -q')G.„(q)nfl(-k, -q, q')G,„(q') , (152) 

k 

where q' = q -|- k The corresponding approximation for the direct correlation 
function is 

AC,,.(k) = Q,-^(k) , (153) 

with Tfei(k) approximated by Eq. (152). 

Several authors [24,27,28,35] have previously obtained a correction to the mean 
field X parameter for a binary homopolymer blend from a Gaussian field theory 
for a nominally incompressible mixture. We now show how the results of these 
earher studies may be recovered from the k = of Eq. (153) for ACjj (k). In a 
binary homopolymer blend, a nonzero value is obtained for f2 -"^ (k) only if all 
of the monomer type indices are equal to a common value i = ii = ■ ■ ■ = i„, in 
which case Q*^"^ is equal to the single-species function VL^^} for the homopolymer 
species a comprised of iV, monomers of type i. To describe such a mixture, we 
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thus adopt a simplified notation in which ilf'' (k) denotes Q^";...j(k). In the 
hmit k ^ 0, we may further simphfy Eqs. (152) and (153) by using the hmiting 
values nf\0,0) = piNf and af^(0,q, -q) = piiV.cjf ^(q, -q), which follow 
from equation (19) for a;, where pi is the number density of homopolymers 
containing monomers of type i. Substituting these limiting forms into Eqs. 
(152) and (153) yields a limit 

lim Aa,(k) ^ ^ / F,(q)G',,(q)G',,(q)F,(q) , (154) 
q 



in which 



^.(q) - ^ (155) 



Eq. (140) may be used to evaluate G in Eq. (154), while setting Jli2(k) = 
^21 (k) = for a homopolymer blend. The corresponding one-loop contribution 
to the difference 

Ax„ = lim[xa(k)-Xo] , (156) 

k— »0 



where Xa(k) is defined by Eq. (143), is given by 



(Fi - F2)V4 



/ [fii +^2- 2Xof^l^^2] 
k 



I 



[Qi + ^2 - 2xo^^i^^2] 



k 



[fll + ^2 - 2X0^^1^^2] 



All integrals in the above may be cut off at a wavenumber /c ~ 7 in order 
to crudely mimic the assumed wavenumber dependence of G(k). Eq. (157) is 
equivalent to the expression obtained for Xa by Wang [35]. It reduces when 
Xo = to the results obtained earlier by Fredrickson, Liu, and Bates [27,28], 
and when Xo 7^ but bi = 62 to the result of de la Cruz, Edwards, and Sanchez 
[24]. 

In all of these earlier applications of the Gaussian field theory, the effective 
X parameter was obtained by examining the macroscopic composition depen- 
dence of a Gaussian contribution to the Helmholtz free energy of a homo- 
geneous mixture, rather than by taking the k — > limit of a wavenumber 
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dependent quantity, as above. The zcro-wavenumber limit of ACij(k) in a 
homopolymer mixture is given by the derivative 

Un>AC«(k) = -:l-^!^ (158) 



in which AA is the difference between the actual Helmholtz free energy of the 
blend and that obtained in a Flory-Huggins approximation. In a homopolymer 
blend, the Gaussian approximation yields a free energy difference 

A A 1 r 

— ~- J ln[(Qi + Q2-2xo^^iQ2)So] (159) 

|k|<7-l 



in which Jli(k) = piu;i(k.). A straightforward differentiation of Eq. (159) with 
respect to molecular number densities yields an expression identical to Eq. 
(154), except for the replacement of ^2 by throughout. This difference arises 
from our use of a renormalized, rather than bare, one-loop approximation to 
obtain Eq. (154). 

Eq. (153) is considerably more general than Eq. (154), or the results of the 
previous work discussed above, insofar as it applies to a wider range of systems, 
including block copolymer melts, and can be used to calculate Xa{^) ioi q 7^ 0. 
This more general result will be needed in order to use the theory presented 
here to describe fluctuations in block copolymer melts near an order disorder 
transition, in which the pre-transitional fluctuations occur near a nonzero 
wavenumber q* rather than near q = 0. 



20.4 Cutoff- Dependence and Renormalization 



The results of the 1-loop approximation for ACij (k) and A^^ depend strongly 
upon the value chosen for the range of interaction 7"^, which acts as a coarse- 
graining length. In a mixture of Gaussian homopolymers with unequal sta- 
tistical segment lengths, bi 7^ 62, the flrst line of Eq. (157) yields a value of 
^Xa ~ 7^ [27,28,35], with a prefactor that vanishes when bi = 62- This corre- 
sponds to a free energy density of order 7^, consistent with the estimate given 
in Eq. (139) for the magnitude of a generic one-loop diagram. The second line 
yields a contribution proportional to Axa that diverges as XoFlt while the 
third line is converges as 7 ^ 00. In the special case of 61 = 62 but Xo 7^ 
considered by de la Cruz et al [24], Eq. (157) thus yields a result Ax^ ~ 7PXo 
that also depends strongly upon 7. 

Because coarse-grained models cannot predict details of atomistic structure, 
they are potentially useful only for describing universal phenomena arising 
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from long-wavelength composition fluctuations. Both the power-counting ar- 
guments of Sec. 19 and the results of the one-loop approximation indicate 
that corrections to the mean field free energy are dominated by contribu- 
tions arising from fluctuations with wavelengths of order the coarse-graining 
length 7"^, and have values that depend strongly upon the value chosen for 
this length. The loop expansion is thus mathematically well-controlled in the 
limit ^ p, but, by itself, not particularly useful, because of the strong 
dependence of all results upon the value chosen for 7. 



To obtain useful predictions for, e.g., the dependence of long- wavelength com- 
position fluctuations upon temperature near a critical point or order-disorder 

transition, the loop expansion must be renormalized, by a process analogous 
to that originally applied to quantum electrodynamics (QED). By "renormal- 
ization" , I mean here a procedure in which contributions to the diagrammatic 
expansion for any quantity that depend strongly upon a cutoff wavenumber 
7 are (if possible) absorbed into a redefinition of phenomenological parame- 
ters, such as the electron mass and charge in QED, or the x parameter and 
statistical segment lengths in a polymer liquid. 



An important first step in this direction was recently taken by Zhen-Gang 
Wang [35]. Wang showed that the above one-loop prediction for Axa can be 
divided into a large "ultraviolet divergent" contribution that depends strongly 
upon 7, but that is independent of N or xA^, and a smaller "ultraviolet con- 
vergent" contribution that is independent of 7, but that depends upon N and 
xN. Only the ultraviolet convergent part, which is insensitive to structure at 
short wavelengths, is sensitive to distance from the spinodal, and develops an 
infrared (i.e., long- wavelength) divergence as the spinodal is approached. Wang 
proposed that the ultraviolet divergent contribution to Axa be absorbed into 
a redefined x parameter. The resulting x parameter, unlike Xo^ corresponds 
to that which would be inferred from the use of the RPA to analyze neutron 
scattering far from any critical point or spinodal. By reexpressing the the- 
ory in terms of this redefined x parameter, Wang was then able to construct 
a renormalized Hartree theory of critical fluctuation whose predictions are 
completely independent of 7. Further refinement of this renormalization pro- 
cedure will be needed for the development of a fully consistent renormalized 
theory of fluctuation effects in block coplymer melts, or of more sophisticated 
theories of critical phenomena in binary blends. This is beyond the scope of 
the present article, but will be pursued elsewhere, using the cluster expansion 
presented here as a starting point analogous to the divergent unrenormalized 
perturbation theory of QED. 
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21 Field Theoretic Approach in Canonical Ensemble 



Most previous applications of the field theoretic approach to polymer liquids 
have started by applying the Edwards transformation to the canonical, rather 
than grand-canonical, partition function [24,25,26,27,28,29,30,31,32,33,34]. The 
choice of ensemble should not effect any exact results or systematic expansions, 
it can affect intermediate steps, approximate results, and the reasoning needed 
to derive diagrammatic rules. 

Application of the Edwards transformation to the canonical partition function 
ZM_[h] of a system containing a fixed number of each species of molecule 
yields a functional integral 

Z[h] j D[J] e^l-^-'^l (160) 

analogous to Eq. (29) for S, in which 

L[J, h]= In Zm-^Ur/jaj . (161) 

k 



Here 



lnZ[^] =^M„ln 



(162) 



is the logarithm of the canonical partition function Z[h] for an ideal gas in a 
field h — h + iJ, and Za[h] is the partition function for a single molecule of 
type a in such a field. 

The most important technical difference between the canonical and grand- 
canonical formulation is that Eq. (162) gives In Z as a sum of terms propor- 
tional to the logarithm Inz^, whereas Eq. (33) gives lnS[/i] = Z^a'^a^al^] as 
a corresponding sum of terms proportional to Za[h] itself. As a result, func- 
tional derivatives of lnZ[/i], which arc needed to construct a diagrammatic 
perturbation theory, are given by sums 



^('^)(l,...,n) = 



1 (^"lnZ[fe] 
V 5h{l) ■ ■ ■ 5h(n) 
Ep„V;i")(l,...,n) 



(163) 
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in which 



n) 



(5" In Za[h] 



(164) 



Sh{l)---Sh{n) 



is an ideal-gas intramolecular cluster function for molecules of species a, which 
is related to the corresponding correlation function uja by a cumulant expan- 
sion. Diagrammatic expansion of the functional integral representation of Z[h] 
thus yields diagrams similar to those obtained here, except for the use of ^ 



In the canonical ensemble formulation in canonical ensemble, all results are 
obtained as explicit functions of molecular number densities, rather than of 
chemical potentials or activities that are known only implicit functions of 
number density. For this reason, the canonical ensemble formulation some- 
times provides the shortest route to useful explicit results at the one-loop or 
Gaussian level. In particular, it was used by Barrat and Fredrickson to obtain 
a very direct derivation of Eq. (150) for the intramolecular correlation function 
a;(^^(q) in a diblock copolymer melt. 

The main disadvantage of the canonical formulation is that the reasoning 
required to develop systematic diagrammatic rules, or to justify various topo- 
logical reductions analogous to those discussed here, is somewhat more com- 
plicated in the canonical formulation. Special rules are required to interpret 
the thermodynamic limit of the Fourier representation of diagrams of ^ ver- 
tices in a homogeneous liquid that are not needed for diagrams of Q vertices, 
in which the expression for the function associated with a vertex depends on 
whether particular subsets of its wavevector arguments have a vanishing vec- 
tor sum. These complications motivated the use of grand-canonical ensemble 
in this article. A full discussion canonical formulation, and of the relationship 
between it and the grand-canonical formulations, will be given elsewhere. 



22 Mayer Cluster Expansion 

Building on the development by Ladanyi and Chandler of a Mayer cluster 
expansions for fluids of rigid molecules [4,8], Chandler and Pratt [5] have 
given Mayer cluster expansions of one- and two-molecule correlation functions 
for interaction site models of fluids of non-rigid molecules. In such diagrams, 
bonds represent factors of the Mayer function 



vertices rather than Q vertices. 



fij{r - r') = e 



Uij{T-r') _ 



(165) 



66 



or its Fourier transform, rather than factors of —Uij itself. In this section, we 
derive the Mayer chister expansion for fluids of flexible molecules by reasoning 
somewhat different than that used by Chandler and coworkers, and outline 
the relationship of the Mayer cluster expansion to the perturbative expansion 
developed in the remainder of this paper. 

22. 1 Rules for Valid Diagrams 

Mayer cluster diagrams for fluids of non-rigid molecules may be constructed 
and interpreted according to rules that are broadly similar to those that apply 
to the perturbative cluster diagrams discussed elsewhere in this paper. For 
concreteness, we will discuss only the rules for the construction of coordinate 
space diagrams, which represent integrals over monomer positions. 

In both Mayer diagrams and perturbative cluster diagrams, each "v vertex" 
with n associated circles represents a function ) , and each h bond 

attached to two circles represents a function 6(r, r'). In both types of diagram, 
black field circles are associated with integration variables, and white circles 
with parameters. In both types of diagram, adjacent circles may not be directly 
connected by more than one bond. 

Mayer cluster diagrams differ from perturbative cluster diagrams in that: 

1) The bonds in Mayer diagrams are / bonds. 

2) Any number of / bonds may be attached to any vertex field circle, and 
bonds may be attached to root as well as field circles. 

3) Each vertex in a Mayer diagram is associated with a particular molecular 
species a, and each vertex circle is associated with a specific physical site 
on a molecule of that species. Different circles on the same vertex must all 
be associated with distinct sites on the corresponding molecule. 

This last rule is more restrictive than the corresponding rule for the construc- 
tion of perturbative cluster diagrams, which require only that each vertex circle 
be associated with a specific "type" of monomer, but not necessarily with a 
specific site on a specific species of molecule, and do no prohibit different sites 
from being associated with the same type. 

To establish a notation for Mayer diagrams similar to what we have used for 
perturbative cluster diagrams, we will herafter restrict each "type" of monomer 
in any Mayer cluster diagram to include only monomers that occupy a specific 
site on a specific species of molecule. For Mayer clusters, we thus abandon the 
more flexible deflnition of type that is useful in the interpretation of perturba- 
tive cluster diagrams. With this convention, each of the monomer type indices 
ii, . . . ,in on an intramolecular correlation function • • • ,i"n) may 
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thus be uniquely associated with a specific site on a molccTilc of type a. With 
this convention, we need not distinguish between the function . . . ,n) 

for a specific species of molecule and the function . . . ,n) that is de- 

fined in Eq. (17) by a sum over species, because the convention implies that 
the sum in Eq. (17) cannot have more than one nonzero term. When the list of 
arguments of an intramolecular correlation function such as i7*^")(l, . . . , n) in- 
cludes two or more monomers with the same type index, and thus refer to the 
same monomer, we take the value of the correlation function, by convention, 
to include a product of 5 functions that constrain the associated monomer po- 
sitions to be equal. For example, if each molecule of type a in a homogeneous 
liquid contain a single monomer of type i, then Q^"^j(ri, r2) = Pa(^(ri — r2). 



22.2 Grand Partition Function 



The Mayer cluster expansion of S for a fiuid of non-rigid molecules may be 
derived by reasoning closely analogous to that normally applied to atomic 
fiuids [3]. The ratio of the grand-partition function S for a fiuid with a potential 
energy of the form given in Eqs. (3-5) to the corresponding partition function 
S of an ideal molecular gas with Uint — 0, may be expressed as an expectation 
value 

J = (e-^"')7.G. (166) 



where (• • •)i,g. denotes an average value evaluated in the ideal gas state. Here, 
S and S are evaluated with the same values of molecular activities and the 
same h fields. The Boltzmann factor e~^^"* may be expressed as a product 

e-^-* = n e-^- =11(1 + f,,) , (167) 



where // and u index all of the monomers present in the system, and 11^,^ is 
the two-body interaction between monomers /i and v. 

By reasoning closely analogous to that discussed by Hansen and MacDonald 
for a liquid of point particles [3], which is discussed in appendix G, we obtain 
a diagrammatic expansion 



Sum of diagrams of Q vertices, 
and one or more / bonds, with no 
root circles 



(168) 
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in which all diagrams are constructed according to the rules OTitlined in Sub- 
section 22.1. A straightforward application of the exponentiation theorem of 
appendix subsection F.l then yields an expansion of ln(S/S) as a sum of the 
subset of diagrams on the r.h.s. of Eq. (168) that are connected. 

The symmetry number associated with a Mayer diagram may be calculated 
according to the rules outlined in appendix C. In the case of Mayer diagram, 
however, because the circles on a vertex must all be associated with specific, 
distinct sites on a molecule, there is no circle symmetry factor arising from 
equivalent permutations of arbitrary labels associated with field circles of the 
same type on the same vertex. The symmetry factor for a Mayer diagram is 
thus given by the order of the symmetry group of permutations of the vertices 
alone. 



22.3 Correlation Functions 



Starting from Eq. (168) for S, one may obtain expressions for correlation 
functions by a process of functional differentiation and topological reduction 
closely analogous to that carried out in Sees. 10-15, with closely analogous 
results. 

Functional differentiation of In 5 with respect to h then yields an expansion 



5(")(l,...,n) 



Sum of connected diagrams of Q 
vertices and / bonds, with n root 
circles labelled 1, . . . , n 



(169) 



that is closely analogous to Eq. (87). Repeating the reasoning of Sec. (11) 
yields an intramolecular correlation function 



fi(")(l,...,n) = < 



Sum of connected diagrams of 
fl field vertices, / bonds, and 
one Cla root vertex with n 
root circles labelled 1, . . . , n 



(170) 



analogous to Eq. (93). In these diagrams, both field and root circles may be 
connected to any number of bonds, but no more than one bond may connect 
any pair of circles. 



The vertices in Eq. (169) and (170) may be renormalized by a procedure closely 
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analogous to that used in Sec. 13. Renormalization of the vertices in Eq. (170) 
yields an expansion 



Sum of connected diagrams 
containing any number of fl 
field vertices and / bonds, 

and a single Qa root vertex } (1^1) 
containing n root circles 
labelled 1, . . . , n, with no 
articulation field vertices 



Q(-)(l,...,n) = < 



which is analogous to Eq. (98). Renormalization of Eq. (169) yields 

5(")(l,...,n) = Q(")(l,...,n) 

Sum of all connected diagrams 
of Q vertices and / bonds 
with n roots circles labelled 
1 , . . . , n on two or more root 
vertices, with no articulation 
vertices 



(172) 



which is analogous to Eq. (100). Equivalent expressions for both the in- 
tramolecular and collective two-point correlation functions have been given 
previously by Chandler and Pratt [5]. In each expansion, the only differences 

between the expansion in diagrams of / bonds and a corresponding expansion 
in diagrams of — [/ bonds are those resulting from the differences in the rules 
for the construction of valid diagrams. 



22.4 Relation of Perturbative and Mayer Expansions 



To relate any diagram of / bonds to a corresponding set of diagrams oi —U 
bonds, it is useful to consider an intermediate type of diagram in which each 
circle on a vertex represents a distinct site on the corresponding molecule, as 
in a Mayer cluster diagram, and in which pairs of circles may be connected by 
any number of — [/ bonds. To relate any diagram of / bonds to a corresponding 
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sum of such diagrams of —U bonds, we Taylor expand each / bonds as 



Applying this to all of the / bonds in a Mayer diagram yields a sum of integrals, 
each of which may be represented as diagram oi —U bonds and Q or Q vertices 
in which any number S^,^ > 1 may connect each pair of adjacent circles /i and 
v. 

When comparing an expansion in diagrams of / bonds to a corresponding 
expansion in diagrams of —U bonds, attention must payed to the conventions 

used to evaluate the interaction [7,^,^ of a monomer with itself. In Eq. (4), a pre- 
factor of 1/2 is included to guarantee that the interaction energy arising from 
each pair of distinct monomers is counted only once. However, this factor also 
implicitly yields a contribution Uii{r — r' — 0)/2 for the self-interaction energy 
of any monomer of type i. To maintain consistency with this convention, we 
must take /^^ — Q-Uu(r=r )/2_]^ ^^^^ each / bond that represents the interaction 
of a monomer with itself. Such bonds are depicted as closed loops. No more 
than one such loop may connect any vertex circle to itself. 

The algebraic expression associated with a diagram of —U bonds and Vl or 
Q vertices that is obtained by Taylor expanding the factors of / in a Mayer 
diagram may be expressed as an integral divided by a combinatorical factor 



where Sy is the vertex symmetry factor for the original Mayer diagram, B^^ 
is the number of —U bonds connecting circles ^ and for any n > i^, and 
= Y^i^Bf^fj, is the total number oi —U bond loops that are connected at 
both ends to the same vertex circle, or the same monomer. The factor of 2'' 
in Eq. (174) is the direct result of our use of a self interaction U^^{v = r')/2, 
which results in this case in a factor of [U^^ijc = r')/2]^ for each vertex circle 
that is connected to itself hy B —U bonds, contributing a factor of 2^ to the 
combinatorical factor in the denominator. 

We may obtain the same intermediate representation from perturbative cluster 
diagrams oi —U bonds by considering perturbative cluster diagrams in which 
(as in Mayer diagrams) each "type" of monomer corresponds to a specific site 
on a specific species of molecule. In general, valid diagrams oi —U bonds may 
contain several circles of the same "type" on a single vertex. In this variant 
of the theory, however, circles of the same type on a single vertex necessarily 
represent the same site on a molecule. To make this equality explicit, we 
may modify our graphical rules in this case by superposing all of circles that 




(173) 



Sv2'- n 



(174) 
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(a) (b) 



Fig. 13. Example of the graphical superposition in a perturbative cluster diagram 
of circles on the same vertex that represent the same site on a molecule. The ex- 
ample shows two equivalent ways of drawing an interaction between two diatomic 
molecules, with sites labelled "1" and "2". In both diagrams, the integer labels 
near each field or root circle indicates the site on the molecule associated with that 
circle. Diagram (a) is a valid perturbative cluster diagram oi —U bonds, in which 
exactly one bond intersects each field site, and none for root sites. Diagram (b) is 
an equivalent diagram in which each circle represents a distinct site, but in which 
multiple bonds may intersect either field or root sites. 

represent the same site on the same vertex. An example of this superposition 
is show in Fig. (13). Superposition of a set of equivalent circles that includes 
one or more root circles produces a root (white) circle, with which we associate 
a fixed position. Superposition of a set of n equivalent field circles produces a 
field (black) circle attached to n —U bonds. 

Each diagram of —U bonds that is obtained by this process of superposing 
equivalent circles is identical to one that may be derived from a base diagram 
of / bonds by Taylor expansion of the / bonds. Furthermore, the symme- 
try number given in Eq. (174) for such a diagram of —U bonds, which was 
obtained above by Taylor expansion of / bonds in a Mayer diagram, may 
be seen to be identical to that obtained in appendix subsection C.3 for the 
corresponding perturbative cluster diagram. Every Mayer cluster diagram of 
/ bonds is thereby shown to be equivalent to a corresponding infinite set of 
perturbative cluster diagrams of —U bonds. 



23 Discussion 

The analysis presented here was initially motivated by my desire for a set 
of systematic rules for the application of diagrammatic methods to Edwards' 
formulation of polymer statistical mechanics. In the process of searching for 
such rules, the various diagrammatic expansions obtained by starting with 
a Wick expansion of the field theory were found to be perturbative variants 
of the Mayer cluster expansions developed by Chandler and co-workers for 
interaction site models of molecular liquids[5], thus clarifying the relationship 
of the two approaches. 

Several diagrammatic resummations have been applied in order to recast the 
renormalized cluster expansion into more convenient forms. The renormalized 
expansion given in Sec. 15 provides a particularly convenient starting point 
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for the analyis of coarse-grained models of polymer liquids, in which collec- 
tive cluster functions are expressed as functionals of the true intramolecular 
correlation functions in an interacting fluid, and of a screened interaction. A 
one loop approximation for this expansion yields results consistent with those 
obtained previously from Gaussian field theory, except for the replacement 
of ideal-gas intramolecular correlations by exact intramolecular correlation 
functions throughout. Specifically, the one-loop calculation given in Sec. 20 
recovers renormalized versions of a result for the two-point intramolecular 
correlation function found previously for diblock copolymer melts by Barrat 
and Predrickson, and for the long- wavelength limit of the effective x parameter 
found in several previous studies of binary homopolymer mixtures. 

One result of this analysis that does not have an obvious antecedent is the 
analysis of a wavenumber dependent collective Ornstein-Zernicke direct corre- 
lation function Cij{k) for a fiuid of non-rigid molecules, and for a corresponding 
wavenumber dependent Flory-Huggins parameter Xa{k) in the incompressible 
limit. The diagrammatic expansion of C(l, 2) obtained here reduces in the ap- 
propriate limit to a well-known expression for the direct correlation function 
of an atomic liquid. The zero-wavenumber limit of a one-loop expansion for 
this quantity has been shown to yield an effective x parameter identical to 
that obtained by other authors by considering the composition dependence 
of the Helmoltz free energy of a homogeneous blend within a Gaussian the- 
ory. The wavenumber dependent expression will be needed in any application 
of the theory to composition fiuctuations in diblock copolymer melts near 
an order-disorder transition, since this will require an accurate description of 
fluctuations with wavenumbers near a nonzero critical wavenumber q*. 

The power-counting arguments and one-loop calculation of Sees. 18 - 20 are 
an attempt to clarify the conditions for the validity of a loop expansion of the 
perturbative cluster expansion for a coarse-grained model of a dense multi- 
component polymer fluid. The loop expansion is found (not surprisingly) to 
provide a rapidly convergent expansion only for very coarse-grained models, 
in which the coarse-grained statistical segment length is much larger than the 
polymer packing length, so that coarse-grained monomers strongly overlap. 

Even in this limit, in which the loop expansion converges, the magnitude of 
the corrections to mean field theory generally depend strongly upon the value 
chosen for a coarse-graining length. As a result of this cutoff dependence, 
it appears that no calculation based on a coarse-grained model can yield a 
meaningful description of the universal effects of long- wavelength composition 
ffuctuations unless it is somehow renormalized, by absorbing cutoff-dependent 
results into changes in the values of a few phenomenological parameters. The 
need for some such form of renormalization has been recognized throughout 
the history of studies of excluded volume affects in polymer solutions, but 
does not seem to have been fully appreciated until recently [35,45] in field 
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theoretic studies of dense polymer liquids. The construction by Wang [35] 
of a renormalized Hartree theory of composition fluctuations in near-critical 
polymer blends has provided an important first step in this direction. A loop 
expansion of the cluster expansion presented here provides a natural basis 
for further work on the development of renormalized theories of the universal 
aspects of composition fluctuations in polymer blends and copolymer melts. 



A Wick Expansion 



In this section, we review the use of Wick's theorem to derive an expansion of 
the functional integral representation of Sanh as an inflnite sum of Feynman 

diagrams. More detailed discussions of diagrammatic perturbation theory arc 
given by Amit [23], Itzykson and Drouffc [44], and Le Bellac [46]. We begin 
by expanding Eq. (69) for Sanh as a Taylor series 




harm • 



The mth term in this series may be expressed as the expectation value of a 
product 

((-f'anh)™) = (-^inh-^inh • • • -^irh)harm (A. 2) 



of numerically identical factors Lanh = -^inh ~ -^anh = ■ ■ ■ Li^^ii , which (for 
bookkeeping purposes) we distinguish by integer labels a = l,...,m. Each 
factor L^^jj may then be expanded as 

4nh= E — / d{aj}L^''-\al,a2,...,ana)J{al)---J{ana) (A.3) 



where a composite index ak is used to identify a fleld J{ak) = Ji^^{Vak) 
that is the /cth factor of J in the expansion of i^aniu corresponding to the kth 
argument of L*^"«\ and where / d{aj} = J d{al) ■ ■ • J d{ana)- 

Substituting Eqs. (A.3) and (A. 2) into expansion (A.l) yields an expression 
for Sanh as an inflnite series in which each term is a multi-dimensional integral 
of the form 

^(fl{^J d{aj}L^''-\al, an,)5J(al) • • • 5 J(an„)) ^ (A.4) 
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where the average (• • •) is calculated using the Gaussian reference distribution 
for J. Each such term may be uniquely identified by a value of m and a list 
of values (ni, 77,2, . . . , n^^), in which the integer identifies which term in 
expansion (A. 3) of the ath factor of Lanh appears in the product within the 
expectation value. 

Using Wick's theorem [23,46,40], we may expand the Gaussian expectation 
value of a product of factors iJ in each term of the form shown in Eq. (A. 4) 
itself as a sum of terms. Each term in this Wick expansion contains a prod- 
uct of factors of K{aj, a'j') — —{5J{aj)5J{a'j') arising from the expectation 
values of pairs of Fourier components, and corresponds to a distinct ways of 
pairing factors of i6J. Each term in the Wick expansion may be represented 
graphically by a labelled diagram of m vertices connected by bonds that rep- 
resent factors of —K. In each such diagram, each L vertex may be labelled by 
a unique integer a — 1, . . . ,m that corresponds to the superscript associated 
with a corresponding factor of L^"|^ in Eq. (A. 2), where m is the total number 
of vertices in the diagram. Each vertex circle on vertex a may also be labelled 
by a label j = 1, . . . ,na, which corresponds to the integer j associated with 
a corresponding factor of iSJ{aj), and with a corresponding argument of the 
function L*^"«)(l, . . . ,71^). Each factor of —{6J{aj)6J{a'j')) in a term in the 
Wick expansion may be represented graphically by a bond between circle j on 
vertex a and circle j' on vertex a'. Wick's theorem guarantees that each term 
in the Wick expansion corresponds to a distinct pairing of factors of i6J{aj), 
so that all of the completely labelled diagrams generated by this expansion 
are topologically distinct. 



B Diagrams and Integrals 

The integral /(F) associated with any diagram F may be intepreted as either a 
coordinate-space integral, defined as an integral with respect to the positions 

associated with all vertex field circles in the diagram, or as a corresponding 
Fourier sum (for a finite system) or integral (in the limit V 0). In any 
representation, /(F) may be obtained by the following prescription: 

(1) For each v vertex with n associated vertex circles with associated monomer 
type indices i, introduce a factor ^^(ri, . . . , r„) into the integrand of a co- 
ordinate space integral, VviCki, . . . , k„) into a corresponding Fourier sum 
for a finite system, or (27r)^5(ki -|- ■ ■ ■ + k„)t;j(ki, . . . , k„) into a Fourier in- 
tegral for an infinite homogeneous system, where ri, . . . , r„ or ki, . . . , k„ 
are sets of positions or wavevectors associated with the vertex circles. 

(2) For each bond, introduce a factor of 6jj(r, r') into the coordinate space 
integral, V"6jj(k, k') into a Fourier space sum for a finite system, or 
(27r)^(5(k + k')6ij(k, — k) into the Fourier integral in an infinite homoge- 
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neous system, where (i, r) and (j, r') or {i, k) and {j, k') are the arguments 
associated with the attached circles. 
(3) For each vertex field circle with associated position r or wavevector k, 
introduce coordinate integral / dr, a Fourier sum J^k, or an integral 

The above rules for Fourier sums and integrals are based on the assumption 
that Vi^,...,i^(ki, . . . , k„) and 6ij(k, k') are defined using convention (11) for 
normalized transforms of functions of n variables. 

The value of the coordinate space integral associated with a diagram is gen- 
erally a function of the coordinates and monomer type index arguments asso- 
ciated with its root circles. If a diagram has no root circles, coordinate space 
and Fourier integrals yield identical values. If the diagram has one or more 
root circles, the value of a Fourier integral associated with a diagram is given 
by the Fourier transform, without the normalizing factor of 1 /V used in in Eq. 
(11), of the value of the corresponding coordinate space integral. The Fourier 
sum associated with a diagram in a large homogeneous system may be shown 
to contain one overall factor of the system volume V for each disconnected 
component of the diagram. 



C Symmetry Numbers 

Here, we consider the determination of the symmetry number S{r) used to 
obtain the value /(r)/S'(r) of an arbitrary diagram F. The diagrams consid- 
ered here may have have arbitrary integer labels associated with some or all of 
their vertices and/or field circles. Though the value of the integral /(F) asso- 
ciated with a diagram F does not depend upon whether or how it is labelled, 
the symmetry number S{r) generally does. The Wick expansion of Sanh dis- 
cussed in appendix A yields completely labelled diagrams with no root circles. 
A completely labelled diagram with m vertices is one in which distinct integer 
labels a = 1, . . . ,m are associated with each of the vertices, and distinct labels 
j — 1, . . . ,n are associated with each of the field circles of each vertex with n 
field circles, as described in appendix A. In completely labelled diagrams with 
root circles, such labels are associated only with field circles, but not with 
root circles. Each root circle is instead uniquely uniquely identified by its as- 
sociation with a specific external position (or wavevector) and monomer type 
argument. An incompletely labelled diagram is one in which some or all of the 
vertices or field circles are left unlabelled, but in which the labels associated 
with any labelled vertices in a diagram with m vertices are distinct integers 
chosen from the set 1, . . . , m, and the labels associated with any labelled field 
circles on a vertex with field circles are distinct integers chosen from the 
set 1, . . . ,na- An unlabelled diagram is one with no such vertex or field circle 
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labels. 



Two completely labelled diagrams are considered equivalent, or topologically 
distinct, if they contain the same number of vertices and free root circles, have 
corresponding numbers of field and root circles on corresponding vertices and 
the same arguments associated with corresponding root circles, and have the 
same pairs of labelled circles connected by bonds. Two incompletely labelled 
diagrams are considered equivalent if it is possible to add missing labels to 
their field circles and/or vertices so as to create equivalent completely labelled 
diagrams. 

The symmetry factor for any completely labelled diagram T is given by 

S(r)^mWn2\---nJ. (C.l) 

where m is the number of vertices in F and Ua is the number of field circles on 
vertex number a. The inverse of S{r) is the overall combinatorical prefactor 
of a term of the form given in Eq. (A. 4), in which the factors of n\ arises 
from expansion of the exponential, and each factors of rijl arises from Taylor 
expansion of one factor of L("-). The same prefactor appears in the value of 
each of the diagrams that is obtained using Wick's theorem to evaluate such 
a term. The Wick expansion thus yields an expression for Sanh as a sum of 
the values of all topologically distinct fully labelled diagrams of L vertices and 
—K bonds with no root circles. 

The value of an incompletely labelled (or unlabelled) diagram F is defined 
to be equal to the sum of the values of all incquivalent completely labelled 
diagrams that may be produced by attaching missing labels to the unlabelled 
vertices and field circles of F. Because the sum of a set of a inequivalcnt 
unlabelled diagrams is thus the same as the sum of a corresponding set of 
inequivalent fully labelled diagrams, it is generally not necessary to distinguish 
between the use of labelled and unlabelled diagrams in expressions that give 
the value of a physical quantity as a sum of the values of a set of all inequivalent 
diagrams that satisfy some restriction. For example, we may describe the Wick 
expansion of In Sanh as a sum of all distinct unlabelled diagrams of L vertices 
and —K bonds with no root circles, or as a corresponding sum of completely 
labelled diagrams. 

C.l Permutation Groups 

To identify symmetry factors for incompletely labelled and unlabelled dia- 
grams, it is useful to introduce a bit of group theory. For the class of diagrams 
considered here, we are interested in groups of permutations of the labels as- 
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sociatcd with the vertices and field circles in a diagram. If A is such a group of 
permutations, let \A\ denote the order of group A. Let the notation Ai = AjA^ 
indicate that A^ is the permutation that is obtained by first applying permu- 
tation Ak and then permutation Aj to the labels of a diagram. To describe 
permutations of labels, it is useful to imagine that we begin with a draw- 
ing of the diagram of interest to which distinct integers have been associated 
with some or all of the vertices and field circles in some specific way, so that 
permutation operations can be described by reference to this initial labelling. 

The group of all allowable permutations of the labels of a completely labelled 
diagram F of m vertices, which we will call -P(r), is the Cartesian product of 
the group of all ml possible permutations of the vertex labels, and of m groups 
of all possible permutations of the labels of field circles around each vertex. 
Note that the symmetry factor for any completely labelled diagram, given in 
Eq. (C.l), is equal to the order of this group: 

S{T) = \P{T)\ . (C.2) 

Let r be an incompletely labelled diagram, and let F' be a valid completely 
labelled diagram that can be obtained by completing the labelling of F. The 
vertex and circle labels that are present in F are the "original" labels of F, 
and those that are added to create F' are "added" labels. Let -P(F) for an 
incompletely unlabelled diagram F denote the group of all allowed permuta- 
tions of labels of the corresponding completely labelled diagram F', so that 
P(F) = -P(F'). Let H{T) denote the group of all possible permutations of only 
the added labels of F', which is a subgroup of -P(F). For a diagram with m 
vertices, 

\H{r)\^m'\n'M\---n'J (C.3) 

where m' is the number of unlabelled vertices in F and n'j is the number of 
unlabelled field circles on vertex number j of F. 

Let G{T) be the group of permutations of the added labels of the completely 
labelled diagram F' that, when applied to F', produce diagrams that are equiv- 
alent to F'. G(F) is referred to as the symmetry group of F. 

The following standard results about subgroups, due to Lagrange [47], will 
be needed in what follows: Let A — {Ai, . . . , A\a\} be a group and B — 
{Bi, . . . , B\B\} be a subgroup of A. A left coset Q of B within A is defined to 
be a subset of A of the form 

Ci = {B,A,B2A,...,B\B\Ai} (C.4) 
which is obtained by multiplying some specific element A^ of A by every el- 
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cmcnt of subgroup B. The number of elements in any left coset of subgroup 
B is equal to (i.e., every element of a left coset may be shown to be dis- 
tinct), and the set of all distinct left cosets of subgroup B form a partition 
of the elements of A into non-intersecting subsets, each containing \B\ ele- 
ments. The same statements also apply also to so-called right cosets, of the 
form {AiBi, . . . , AiB\B\}- If there are N distinct left or right cosets, it follows 
that \A\ = A^|-B|, and thus that the order of any subgroup of A must be 
an integer divisor of \A\. The integer = 1/11/ 15| is known as the index of 
subgroup B within A. 



C.2 Unlabelled and Incompletely Labelled Diagrams 



The value of any incompletely labelled diagram F may be expressed as a ratio 



where A^(r) is the number of inequivalent ways of completing the labelling of 
r. The corresponding symmetry number is thus 

s(r) = ^ . (C.6) 



for any incompletely labelled diagram. 

Let r be an incompletely labelled diagram, and let V be a specific completely 
labelled diagram that may be obtained by completing the labelling of F. The 
set of all possible ways of completing the labelling of F may be generated by 
applying to F' every element of the group (F) of all possible permutations of 
the added labels. The resulting set of fully labelled diagrams (or corresponding 
permutations) may be divided into subsets such that the members of each 
subset are related by elements of the symmetry group G{T), and are thus 
equivalent, while diagrams in different subsets are all inequivalent. 

All the elements of such a set of equivalent diagrams may be obtained by start- 
ing with any one member of the subset, which may be obtained by applying 
some permutation Hj to the labels of F', and then applying every element of 
G{r) to that diagram. If the elements of G(F) are denoted Gi, . . . , G\g\: any 
such set of equivalent diagrams may thus be obtained by applying the set of 
permutations 

{GiH„G2H„...,G\G\Hj} (C.7) 
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to r' for some Hj. Each set of equivalent permutations of the labels of V thus 
corresponds to one of the left cosets of G{T) within H(V). This implies that 
there are |G(r)| equivalent diagrams in each such subset. It also follows that 
the number of such subsets, which is equal to the number of inequivalent ways 
of completing the labelling of F, is given by the index of symmetry group G'(r) 
within H{V'): 



Combining this with Eq. (C.6) yields a symmetry factor 

S{T)^\G{T)\L{T) (C.9) 



where 



^ ' ~ \E{^)\ m'\n'^\n'2\---n'J ^ ' ' 



is the index H{T) in -P(r). The factor L{r) is the number of valid ways of 
choosing a set of labels for the labelled vertices and circles of F from among the 
integer labels of F'. If F is completely labelled, L{T) = |-P(F)| and |G(F)| = 1, 
so that we recover our starting point S{T) = |P(F)|. If F is unlabelled, L{T) ~ 
1, and Eq. (C.9) thus reduces to 

5(F) = |G(F)| . (C.ll) 



The symmetry factor of an unlabelled diagram is thus equal to the order of 
its symmetry group. The difference between Mayer diagrams and perturbative 
cluster diagrams (or Feynman diagrams) is the required definition of the group 
of allowed permutations, which includes only permutations of vertex labels in 
Mayer diagrams, but must include permutations of the circle (or bond end) 
labels in perturbative cluster diagrams or Feynman diagrams [48,49]. 



C.3 Vertex and Circle Symmetry 



The calculation of |G(F)| for an incompletely labelled perturbative cluster 
diagram F may be simplified by expressing |G(F)| as a product of factors 
associated with permutations of the vertex labels and with permutations of 
the circle labels. Let Gc{^) be the subgroup G(F) that leave all vertex labels 
unchanged, and involve only permutations of added field circle labels. This 
will be referred to as the circle symmetry group. If we partition G(F) into 
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left cosets of Gc(r) in G{T), the elements of each coset will be related by 
circle permutations alone, and thus correspond to the same labelling of the 
vertices, while permutations in different cosets will correspond to different 
ways of labelling the vertices. 

Let r" be a diagram that is obtained by labelling any unlabelled vertices of 
r and removing the labels from any labelled field circles of F. Let Gv{^), the 
vertex symmetry group of F, be the group of permutations of the vertex labels 
added to F to produce F", so as to produce diagrams of labelled vertices and 
unlabelled circles that are equivalent to F". Each diagram that is generated 
by applying one of the elements of Gy(F) to F" corresponds to a labelling of 
vertices used in one of the cosets of 5'c(F) within S'(F). The order |G'v(F)| 
is thus equal to the index of the circle symmetry group Gc(r) within G{T), 
implying that 

|G(F)H|Gc(F)||Gv(F)| . (C.12) 

The circle symmetry group Gc(r) of an incompletely labelled diagram may 
be generated from a few types of operations, all of which involve pairs of ver- 
tices that are connected by more than one bond, or vertices that are directly 
connected to themselves by one or more loops: Specifically, elements of Gc(r) 
may be generated from [48]: i) Permutations of sets of two or more indistin- 
guishable bonds that all connect sites of type i on one vertex to sites of a single 
type j on a second vertex, by exchanging pairs of labels for the corresponding 
pairs of circles, ii) Permutations of pairs of circle labels associated with bonds 
that form loops connecting two labelled field circles on the same vertex, and 
iii) Exchanges of the labels associated with the two ends of bonds that connect 
two labelled field circles of the same type on the same vertex. 

To calculate | Gc (F) | , we thus include a factor of Bai^pj ! for each pair of vertices 
a and /3 that are connected by Bai,i3j bonds between labelled field circles of 

types i and j, including the case a = /? of a vertex a that is connected to 
itself by B^i^aj such loops. We also include an additional factor of 2 for each 
bond that connects pairs of labelled field circles of the same type on the same 
vertex. This gives [48,49] 

|Gc(r)|=2'- n (C.13) 

ai,l3j 

where r = J2ai Bai,ai is the total number of bonds that are connected to two 
labelled field circles of the same type on the same vertex. The product in Eq. 
(C.13) is taken over all distinct types of bonds, by taking only a> (3, with all 
values of i and j for a > /3 and all i> j ior a — (5. 

The Mayer cluster diagram discussed in Sec. 22 can have no more than than 
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one bond between circles of types i and j on any given pair of vertices, where 
each monomer type correspond in this context to a specific site on a particular 
species of molecules. As a result, -Bm,/3j < 1 for all ai and Pj. and |Gc;(r)| = 1. 
For Mayer clusters, the calculation of a symmetry number thus requires only 
the determination of |Gy(r)|. 



D Diagrams with Generic Field Circles 

The definitions of the integral and symmetry number associated with a di- 
agram in the preceding two appendices required that each field circle in a 
diagram be associated with a specified monomer type index. The choice of a 
list of monomer types index values for the field circles was thus taken to be 
part of what defines a particular diagram. With this definition, expansions of 
physical quantities such as correlation functions as sums of infinite sets of dia- 
grams must generally include sets of topologically similar diagrams that differ 
only in the values chosen for these monomer type indices, and that would 
become indistinguishable if these monomer type indices were erased. 

In this section, we define diagrams with generic field circles, in which monomer 
types are not specified for any field circles. Each such diagram represents a 
sum of the values of a set of topologically similar diagrams in which monomer 
type indices are specified for all field circles. Diagrams with generic field circles 
may be either completely labelled or unlabelled. A completely labelled diagram 
with generic field circles is one in which arbitrary integer labels are assigned 
to all of the vertices and field circles, with vertices labelled a = 1, . . . ,m and 
field circles around vertex a labelled k = 1, . . . ,na, but in which no monomer 
type is specified for any field circle. The value of a completely labelled diagram 
with generic field circles is defined to be equal to the sum of the values of all 
inequivalent completely labelled diagrams that can be obtained by assigning 
a specific monomer type to each of the field circles. The fully labelled dia- 
grams that are obtained by this process from the same diagram with generic 
field monomers are all inequivalent, because each associates a different fist of 
monomer types with a list of labelled field circles, but they share a common 
symmetry factor, given by Eq. (C.l). The value of a fully labelled diagram F' 
with generic field circles may thus be expressed ratio J(r')/S'(r'), in 

which /(F') is the unrestricted sum of the integrals associated with diagrams 
obtained by assigning all possible monomer types to each field circle (without 
regard to whether the resulting diagrams are inequivalent, since they neces- 
sarily are), and >S'(F') is the symmetry factor common to all of these diagrams. 
Two completely labelled diagrams with generic field circles are equivalent if 
they contain the same set of connections between pairs of labelled generic field 
circles, in which each field circle is identified by a composite index ak. 
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The value of an unlabelled diagram T with generic field circles is defined to 
be the sum of the values of all fully labelled diagrams with generic field circles 
that can be obtained adding arbitrary labels to all of the vertices and field 
circles in F, while leaving the monomer types unspecified. Two unlabelled 
diagrams with generic field circles are equivalent if it is possible to add vertex 
and field site labels to the two diagrams in such a way as to generate equivalent 
fully labelled diagrams with generic field circles. By repeating the reasoning 
applied previously to diagrams with specified monomer types, we find that 
the value of an unlabelled diagram F with generic field circles is given by a 
ratio /(F)/5'(F), in which 7(F) is the integral associated with a corresponding 
completely labelled diagram F', as defined above, and S'(F) is the order of the 
symmetry group of permutations of the labels of F' so as to produce diagrams 
equivalent to F'. The symmetry group for a diagram with generic field sites 
may be determined by simply treating all field sites as if they shared a common 
"generic" monomer type. 

As an example, consider the symmetry number for the ring diagrams discussed 
in subsection 9.2. Let Rn be the unlabelled ring diagram with n vertices and 
2n generic field circles. Consider a reference labelling of the vertices in which 
successive vertices are labelled 1, 2, . . . , n clockwise around the ring. The n — 1 
diagram has a circle symmetry number Sc = '2 and a vertex symmetry number 
Sv = 1, giving an overall symmetry number 2. The n = 2 diagram has Sc = 2 
and Sv = 2, giving 5*2 = 4. All such diagrams with n > 2 have Sc = 1 and 
Sy — 2n, because the connections between n vertices are invariant under n 
possible cyclic permutations (1— >2, ...n— and under n permutations 
that are obtained by pairwise exchange of vertices i and n+l—i (i.e., reflection 
through a plane that cuts the ring in half) followed by any of the n possible 
cyclic permutations. For any n> 1, we thus find a symmetry number Sn = 2n. 
The value of the integral associated with a ring diagram. The summation over 
all possible values of the 2n monomer type indices needed to calculate I{Rn) is 
accomplished automatically by the matrix multiplication and trace operations 
used in Eq. (81). 



E Derived Diagrams 

In appendix F we consider several procedures in which a set of diagrams is 
graphically derived from some base diagram by applying each of a specified set 
of graphical modifications to the base diagram. In appendix subsection F.2, 
we consider diagrams that can be obtained by adding one root circle to any 
vertex of a base diagram, or by inserting a vertex with one root circle into 
a bond. In subsections F.3 and F.4 we consider set of diagrams that may be 
obtained by replacing specific sets of vertices or bonds of a base diagram with 
any member of some specified set of subdiagrams. In this appendix, we derive 
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a lemma that is useful in the evaluation of sums of such derived diagrams. 

Let r be an unlabelled base diagram, and let V be a completely labelled base 
diagram that is obtained by completing the labeUing of F. Let A be the set 
of diagrams that may be obtained by applying any of a specified set B of 
graphical modifications to V. Different elements of B will be distinguished by 
an integer index a, (3 = 1, . . . , |B|. Let be the labelled derived diagram that 
is obtained by applying modification a to T'. Let A^, be the unlabelled derived 
diagram that is obtained by removing all vertex and field circle labels from 
A'^. For example, if A is the set of diagrams that may be obtained by adding 
a root circle with specified spatial and monomer type arguments to any vertex 
of F', and a is taken to be the label in F' of the specific vertex to which a root 
circle is added, then A^ is the labelled diagram that is obtained by adding 
one such root circle to vertex a of F'. If a graphical modification a adds field 
vertices or field sites to F' (as when a new vertex is inserted into a bond) the 
added vertices and sites are left unlabelled in A^. Let G*{Aa) be the group 
of permutations of the labels of A'^ (which are in one-to-one correspondence 
with those in F') so as to create diagrams that are equivalent to A^ 

We assume that set of B of allowed graphical modifications are defined such 
that: 

i) The symmetry group G*(Aq) is a subgroup of the symmetry group 5'(F) of 
the base diagram. 

ii) Any diagram that can be obtained by applying an element of the symmetry 
group G(F) to A'^ is equivalent to one that may be obtained by applying 
one of the elements of B to F'. 

Assumption (i) requires, in essence, that the allowed graphical modifications 
not introduce any symmetries that were not already present in the base dia- 
gram, though it may reduce the symmetry of a diagram. This is satisfied for 
the example of the addition of a root circle to any one vertex, for which the ad- 
dition of a root circle to a vertex makes that vertex distinguishable from every 
other in the diagram, thus potentially lowering the symmetry of the diagram. 
This example may also be shown to satisfy assumption (ii). Assumption (ii) 
would be violated, however, by a set of modifications that allowed a root site 
to be added only to some arbitrarily chosen subset of "modifiable" vertices, 
identified by their labels in F', if this subset is not invariant under application 
of the elements of symmetry group G'(F). Assumption (ii) will generally be sat- 
isfied whenever the specified set of graphical modifications can be described 
without reference to an arbitrary choice of how the diagram F' is labelled. 

Theorem: Let N{Aa) be the number of distinguishable completely labelled 
diagrams that can be obtained by applying elements of B to a labelled base 
diagram F', but that can be reduced to the same unlabelled diagram A^, by 
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removing all labels. Given the assumptions stated above, 



i.e., N{Aa) is given by the index of G*{Aa) within G{r). 

Proof. Let C(Aa) be the set of |G'(r)| diagrams that arc obtained by applying 
all of the permutations in G{T) to A^. This procedure generates all possible 
ways of labelling so as to create diagrams that are labelled in a manner 
consistent with that used in F', but that (by construction) can be reduced to 
Aq by removing all of the added labels. Assumption (ii) states that all of the 
diagrams in this set can be generated by applying graphical modifications in 
set B to r'. To determine the number A^(r) of distinguishable diagrams in 
C(Aa), we divide C(Aq) into subsets, such that diagrams in the same subset 
are equivalent, and must thus must be related by elements of the symmetry 
group G{Aa)j while diagrams in different subsets are inequivalent. The number 
of such subsets, and thus the number of distinguishable diagrams, is given by 
the index of G*{Aa) in G'(r). 



F Fundamental Lemmas 



In this section, we outline or cite proofs of several lemmas about diagrams are 
needed to rigorously justify the topological reductions discussed in the main 
text. Each of these is a generalization of one of those originally presented by 
Morita and Hiroike [1], and reviewed by Hansen and MacDonald [3], for the 
case of point particles. In their original forms, these lemmas applied to Mayer 
diagrams for fluids of point particles, in which each vertex represents a function 
of one position, which is traditionally represented in liquid state theory by a 
circle, and in which no more than one bond may connect any pair of adjacent 
vertices. This appendix presents the generalization necessary for perturbative 
cluster diagrams for molecular liquids, in which each vertex may represent a 
function of several variables, which may itself be a functional of a set fields 
conjugate to the monomer concentrations, and in which multiple bonds may 
connect indistinguishable pairs of circles on adjacent vertices. The counting 
of symmetry numbers relies on the generalized definition of symmetry groups 
and symmetry numbers discussed in appendix C. 
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F. 1 Exponentiation 



The first lemma corresponds to Lemma 3 of Morita and Hiroike and to a 
special case of Lemma 1 of Hansen and MacDonald. It has been called both 
the exponentiation theorem [3] and the linked cluster theorem [40]. 

Theorem 1: Let A be a set of N distinct, connected unlabelled diagrams 
Fi, . . . , Fjv (where N may be infinite), and let 



m 



Li + • • • + r 



N 



fF.ll 



be the sum of values of the diagrams in A. Then 



1+ <^ 



Sum of diagrams comprised of 
components that belong to set 
A, each of which may appear 
repeated any number of times 



> . 



(F.2) 



Proof: The proof of this theorem is identical to that given by Hansen and 
MacDonald [3] (see page 527) for Mayer cluster diagrams or by Huang [40] for 
the Feynman diagram expansion of the free energy (or the generating function 
for connected Green's functions) in a statistical field theory. In this case, the 
required counting of symmetry factors for diagrams with multiple components 
is independent of the nature of the component diagrams. 

Comment: Lemma 1 of Hansen and MacDonald actually contains both Lemma 
4 of Morita and Hiroike, which is the more general form of the theorem, and 
their Lemma 3, which is a special case essentially identical to the theorem 
given above. The more general form of the exponentiation theorem relates the 
exponential of a sum m of a set A of "star-irreducible" diagrams [3] to the 
sum of all diagrams in set A and of all "star products" of diagrams in A. This 
more general result will be referred to here as the star product theorem. 

The star product theorem is needed in the derivation of a density expansion 
for the chemical potential for an atomic liquid (see pages 90 and 91 of Hansen 
and MacDonald). The expansion of the chemical potential is in turn used to 
obtain a density expansion of the Helmholtz free energy. Unfortunately, the 
star product thereom cannot be trivially generalized to the types of diagrams 
used here to describe molecular liquids, in which the vertices represent func- 
tions of more then one monomer position. For this reason, no explicit density 
expansions are given for the chemical potential and Helmholtz free energy in 
this article. Several nontrivial generalizations of the star product theorem will 
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be presented in a subsequent paper, where they wiU be used to derive expan- 
sions of these and several other quantities as exphcit functions of molecular 
number densities. 



F.2 Functional Differentiation 



In this subsection, we derive diagrammatic expansions for functional deriva- 
tives of the values of diagrams in which the functions associated with the 
vertices are themselves functionals of a field hi{r). We consider a general class 
of diagrams in which different functions may be associated with root vertices 
and field vertices, in which all root vertices are v vertices and all field vertices 
are u vertices. We assume that the function associated with a root vertex 
with n field and root circles be an n-th functional derivative 

... n)^ (F3) 



of some generating functional v[h], so that the v-functions obey a recursion 
relation of the form 



as in Eq. (85). In addition, we assume that functional differentiation of the 
function m*^") associated with a field vertex also yields 



so that functional differentiation of a field u vertex yields a root v vertex. 

Below, we consider separately a case in which the factor associated with the 
bonds does not depend on /i, as in an expansion in —U bonds, and a special 
case in which the bonds represent screened interactions that do depend upon 
/i in a specific manner. 



F.2.1 Constant interaction bonds 

Lemma 2. a: Let T be an unlabelled diagram diagram containing any number 
of u field vertices, v root vertices and h bonds, in which there are (n — 1) root 
sites labelled 1, . . . , n — 1. Let the functions u and v be functionals of a multi- 
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component field h that obey Eqs. (F.4) and (F.5), and let bij be independent 
of h. Then 



sr 

Sh{n) 



Sum of unlabelled diagrams that 
can be derived from F by adding 
one root circle labelled n to any 
vertex of F 



> . 



(F.6) 



The diagrams that are obtained by adding a root circle labelled n to one of 
the vertices of F are also interpreted as diagrams of u field vertices, v root 
vertices, and b bonds. 

This theorem is closely analogous to Lemma 2 of Hansen and MacDonald or 
Lemma 1 of Morita and Hiroike. 

Proof. Let F' be completely labelled diagram obtained by adding labels in a 
specific way to all vertices and field circles of F. The desired derivative of F is 
given by 



sr 



1 si{r') 



Sh{n) S{r) Sh{n) 



(F.7) 



where liT') denotes the integral associated with either F or F'. The integrand 
of /(F') is a product of factors arising from the vertices, which are functionals 
of the field, and from the bonds, which (in this case) are assumed to be inde- 
pendent of the field. Use of the chain rule to evaluate the functional derivative 
of this integrand will thus generate a sum in which each term is generated by 
differentiation of a different labelled vertex of F', each of which is associated 
with a diagram in which a root circle is added to that vertex. The derivative 
ST /Sh{n) of a diagram with m vertices is thus given by a sum 



ST 



Sh{n) S{V) i 



:f.8) 



where is the diagram that is obtained by adding one root circle labelled n 
to vertex a. of F'. 

Labelled diagrams and ISl^ that are obtained by differentiating different 
vertices a and /3 of F' are always distinguishable, but may correspond to 
equivalent unlabelled diagrams. Eq. (F.8) may thus be rewritten as 



ST 



1 



Sh{n) 5(F) 



EiV(Aa)/(A:.) 



(F.9) 
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Here, denotes a sum over all values of a that yield inequivalcnt unlabelled 
diagrams, and A^(Aq) is the number of ways of choosing a vertex to which to 
add a root circle so as to obtain a diagram that, upon removal of all labels, 
reduces to one equivalent Aq. Using the Lemma of appendix E, we find that 
N{Ka) = \G{V)\/\G{Ka)l thus yielding the required form 



5h(n) V |G(A„) 



in which |G(Aa)| is the symmetry factor for diagram A^. 



F.2.2 Screened Interaction Bonds 

We now consider a situation in which the factors associated with the bonds 
represent a screened interaction that is also a functional of the field. To ab- 
stract the essential features of a diagram of —G bonds, we consider diagrams 
of —6 bonds for which 

ri(l,2)=«(2)(i^2) + c(l,2) , (F.ll) 



where c(l, 2) is a function that is independent of h. 

Lemma 2. b: Let F be an unlabelled diagram diagram containing any number 
of u field vertices, v root vertices and —b bonds, in which there are (n — 1) 
root sites labelled 1, . . . ,n — 1. Let the functions u and v be functionals of a 
multicomponent field h that obey Eqs. (F.4) and (F.5), and let 6 be a function 
of the form given in Eq. (F.ll), in which c is independent of h. Then 



sr 

6h{n) 



Sum of unlabelled diagrams that can 
be derived from F by either adding 
one root circle labelled n to any vertex 
of F, or by inserting a v vertex with 
one root circle labelled n and two 
field circles into any bond of F 



(F.12) 



Proof. In this case, the required derivative is a sum of terms arising from 

differentiation of different vertices and terms arising from differentiation of 
different bonds. The treatment of terms arising from differentiation of vertices 
is the same as in the case of constant b bonds. To complete the proof, we must 
consider the terms that arise from differentiation of bonds. 
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Let the bonds in the base diagram be indexed by an integer j3, and let A'^ be 
the labelled derived diagram that is obtained by differentiating bond f3. The 
sum of terms arising from differentiation of all possible bonds may then be 
expressed as a sum 



in which denotes a sum over all bonds, or as a sum 

-^X:^(A,)7(A^) (F.14) 



in which J2'f3 represents a sum over choices of bonds that lead to inequivalent 
unlabelled derived diagrams, and -/V(A^) is the number of ways to choose 
a bond so as to obtain unlabelled diagrams equivalent to A^. By repeat- 
ing the reasoning outlined in appendix section E, we find that A'"(A^) = 
|G'(A^)|/|G*(A^)|, and thus that the contribution to 5V/5h{n) of terms arising 
from differentiation of the bonds is of the required form 



thus completing the proof. 



F.3 Vertex Decoration 



In this subsection, we prove a generalization of Lemma 5 of Morita and Hiroike 
or Lemma 4 of Hansen and MacDonald. 

Lemma 3: Let F be an unlabelled diagram containing any number of field u 
vertices, root v vertices and h bonds. Let F contain a specified set of m "target 
vertices" , that each contain exactly n circles with a specified (unordered) set 
of monomer type indices i — {ii, . . . , The target vertices must either be 
all u vertices or all v vertices. The n circles on a root target vertex may be any 
combination of field and/or root circles. Let the membership of the specified 
set of target vertices be unchanged by any isomorphism of F, i.e., by the 
application of any element of the symmetry group of F. 

Let A be a set of inequivalent unlabelled diagrams, each of which contains 
exactly one root vertex that is of species a and contains n root circles labelled 
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1, . . . , n with monomer type indices corresponding to those of the target ver- 
tices. Let 

m(l, . . . , n) = I Sum of diagrams in A | (F.16) 



Let r* be the unlabelled diagram that is obtained by replacing the function 
v(ki, . . . , kn) associated with each target vertex in F by m(l, . . . , n), thereby 
changing the target vertices from u- or v vertices to m vertices. Let B be the 
set of all incquivalcnt unlabelled diagrams that can be obtained by decorating 
every target vertex of F with any one of the diagrams belonging to A, in the 
process shown in Fig. 7, by superimposing the root vertex of a diagram in A 
on each target vertex and replacing the n circles of the target vertex by the n 
root circles of the corresponding diagram in A. If each diagram in B can be 
uniquely reduced to F by removing the pendant diagrams, then 

r* = < Sum of diagrams in B 



Proof: Let F' be a diagram that is obtained by attaching labels to all of the 
vertices and field circles of F in some arbitrary way. Let 

A' = A'(7i,...,7m) (F.18) 



be a diagram that is obtained by attaching a specified list of unlabelled pen- 
dant diagrams {71, . . . , 7^} belonging to A to the m target vertices of F, where 
the value of 7^ specifies the diagram in A that is used to decorate target vertex 
i. In the language of appendix E, each allowed graphical modifications of F' 
corresponds to a different ordered list (71, . . . ,7m) of pendant subdiagrams. 
Let A = A(7i, . . . , 7m) be the unlabelled diagram that is obtained by removing 
all vertex and field circle labels from A'. 

The value of F* may be expressed as a sum 

T.*_ 1 /(AX7i,---,7m)) ^p^g^ 



in which each of the indices 71,..., 7^ may run independently over values 
corresponding to all of the diagrams in A. Each of the resulting lists of values 
for (71, ... , 7m) leads to a distinct labelled derived diagram A' (71, ... , 7m), but 
some lists may may lead to labelled diagrams that correspond to equivalent 
unlabelled diagrams. We may thus rewrite Eq. (F.19) as 

-^(r) 7i,-,7m '^(7l)---'^(7m) 
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where I]^!,...,^^ denotes a sum over all lists of pendant diagrams (71, . . . ,7^) 
that lead to inequivalent unlabelled diagrams A(7i, . . . ,7f„), and A^(A) is the 
number of ways of choosing a list (71, . . . ,7m) so as to produce labelled de- 
rived diagrams A' that reduce to equivalent unlabelled diagrams upon removal 
of all labels. By applying the Lemma of appendix E (which relies upon our 
explicit assumption that the set of target vertices is invariant under the vertex 
symmetry group of F) we find that, for an unlabelled derived diagram A, 



where G*{A) is the group of permutations of the labels of the labels of labelled 

derived diagram A' so as to produce equivalent diagrams. Note that in a 
labelled derived diagram A' with decorated vertices, only the vertices and 
fields circles present in F are labelled, while any field vertices and field circles 
of the pendant subdiagrams are unlabelled. 

Let A" be a diagram that is obtained by completing the labelling of A', by 
adding labels to all of the field circles and field vertices of A'. The symmetry 
group of the unlabelled derived diagram A with pendant subdiagrams on all of 
its target vertices is the Cartesian product of the group G*{A) of permutations 
of the labels of A' and of the internal symmetry groups of pendant subdiagrams 
7i, ■ ■ ■ , 7m- The overall symmetry number <S'(A) is thus given by a product 

m 

|G(A)H|G*(A)iniG(7.)l ■ (F-22) 

i=l 



By combining Eqs. (F.20), (F.21), and (F.22) we obtain an expression for F* 
as a sum 

P* ^ ^ /(A(7i,...,7m)) . . 



of the values of all of the distinguishable diagrams in B. 



F.4 Bond Decoration 



This theorem is analogous to Lemma 6 of Morita and Hiroike or Lemma 5 of 
Hansen and MacDonald. 

Lemma J^: Let F be an unlabelled diagram containing v vertices and h bonds. 
Let there be exactly m b bonds in F that connect circles of monomer species i 
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and J, which will be referred to as target bonds. Let A be a set of inequivalent 
unlabelled diagrams that each contain exactly two free root circles, labelled 1 
and 2, with monomer types i and j, respectively. Let 



Let r* be the unlabelled diagram that is obtained by replacing the function 
b associated with each target bond in F by m, thereby changing the target 
bonds from b bonds to m bonds. Let B be the set of all inequivalent unlabelled 
diagrams that can be obtained by replacing each target bond F with a diagram 
belonging to A, by superimposing the free root circles of each diagram in set 
A onto the circles that terminate the target b bond. Then 



Proof: The proof is essentially identical to that given for the case of ver- 
tex expansion, except that the allowed set of graphical modifications involve 
replacement of a bond, rather than a vertex, by any of a specified set of sub- 
diagrams. Note that the simple analogy between the proofs of the vertex and 
bond decoration theorems relics upon our use of the definition of symmetry 
groups introduced in appendix C, in which the allowed symmetry operations 
are taken to include permutations of field circle labels as well as permutations 
of vertex labels. 



G Mayer Cluster Expansion for Moleculeir Liquids 

In this appendix, we outline the derivation of a Mayer cluster expansion of S 
for a fluid of permanent (i.e., nonreactive) but non-rigid molecules by reasoning 
closely analogous to that usually used to derive an activity expansion of S 
for an atomic liquid. The analysis starts from Eq. (167) for the ratio S/S. 
As in the expansion of S for an atomic liquid, we may expand n(l + fap) 
into an inflnite set of terms, each of which consists of one or more factors 
of fai3 for different pairs of monomers. Each term in this expansion may be 
represented schematically by a diagram in which monomers are represented by 
black circles, an / bond is drawn between each pair of circles a and f3 for which 
there exists a factor of fap in the term of interest. Circles that correspond to 
monomers that are part of the same molecule are drawn as fleld circles around 
a corresponding vertex. In such diagrams, each circle must be attached to at 
least one bond. 

The identities of the monomers in such a diagram may be specified by labelling 




(F.24) 




(F.25) 
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each vertex with a species index a and a molecule index m = 1, . . . , Ma, and 
each field circle by a site index. The resulting diagrams will be referred to here 
as a Mayer diagrams with labelled molecules. Two diagrams with labelled 
molecules are equivalent if they correspond to the same set of connections 
between specific pairs monomers, and thus to the same product of /-functions. 
The value of each such diagram is defined to be the expectation value of the 
corresponding product of /-functions in an ideal-gas with fluctuating numbers 
of molecules and a specified set of molecular activities. To describe systems 
with fiuctuating numbers of molecules, molecules of species a are labelled 1 
to Ma in microstates in which the system system contains exactly Ma such 
molecules. When evaluating expectation values, the value of a product of /- 
functions involving monomers on a specified set of molecules is taken to vanish 
in microstates in which any of the specified molecules is not present in the 
system. The ratio S/S is then equal to the sum of all topologically distinct 
Mayer diagrams with labelled molecules. 

For each Mayer diagram V with labelled molecules, we also define an unla- 
belled Mayer diagram F that is obtained by removing the molecule label m 
from each vertex in F', while leaving the molecular species index a for each 
vertex and the site index s for each vertex field circle. The value of F is defined 
to be the sum of the values of all topologically distinct Mayer diagrams with 
labelled molecules that may be obtained by assigning positive integer labels to 
all of the vertices in F, with distinct labels for all vertices of the same species. 
By construction, the ratio S/S may thus also be expressed as a sum of all 
topologically distinct unlabelled Mayer diagrams. 

The value of an unlabelled diagram F may be expressed as a ratio /(F)/S'(F), 
in which J(F) is the sum of values of the infinite set of all Mayer diagrams 
with labelled molecules that may be obtained by associating molecular labels 
with each of the vertices in F, subject only to the requirement that all vertices 
of the same species have distinct labels, but without requiring that different 
diagrams in the set all be topologically distinct. The symmetry number S(T) 
is given by the order of the group of permutations of the molecule labels of F' 
that generate equivalent diagrams. 

The value of the integral /(F) may be expressed as an integral with respect to 
a set of monomer positions whose integrand is a function of a set of /-functions 
times a probability of finding monomers of specified types at a specified set of 
locations, without any restriction on the labels of the molecules to which the 
monomers are attached except for the requirement that different vertices of the 
same species in F represent distinct molecules. In grand-canonical ensemble, 
the positions of monomers on distinct molecules within an molecular ideal 
gas are statistically independent, and so the probability distribution required 
in the integrand of /(F) may be expressed as a product over single-molecule 
probability densities. The probability of finding a set i of monomers of a 
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molecule of type a in a set of positions r in an ideal gas reference state is 
given by the function ^^^"^(r). The integrand of /(F) may thus be obtained by 
simply associating a factor of Cl with each vertex and a factor of / with each 
bond in F. 
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